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1 Introduction 

One of the most challenging problems in particle physics is to understand hadron dy- 
namics and spectroscopy in terms of the confined quark and gluon quanta of quantum chro- 
modynamics, the fundamental theory of the strong interactions. A central goal is to compute 
detailed hadronic properties, such as moments, structure functions, distribution amplitudes, 
transversity distributions, elastic and transition form factors, and the excitation dynamics of 
hadron resonances from first principles; i.e., directly from the QCD Lagrangian. The most 
successful theoretical approach thus far has been to quantize QCD on discrete lattices in 
Euclidean space-time [1]. Lattice numerical results follow from the computation of frame- 
dependent moments of distributions in Euclidean space; however, dynamical observables in 
Minkowski space-time, such as the time-like hadronic form factors, are not obtained directly 
from Euclidean-space lattice computations. Dyson- Schwinger methods have led to many im- 
portant insights, such as the infrared fixed-point behavior of the strong coupling constant [2]; 
however, in practice, these analyses are limited to ladder approximation in Landau gauge. 

In principle, one could calculate hadronic spectroscopy and wavefunctions by solving for 
the eigenvalues and eigenf unctions of the QCD Hamiltonian: -f^|\E') = -El^P) at fixed time t. 
However, this traditional method - called the "instant form" by Dirac [3], is plagued by com- 
plex vacuum and relativistic effects. In contrast, quantization at fixed light-front (LP) time 
T = t + z/c - the "front-form" of Dirac [3] - provides a powerful boost-invariant nonpertur- 
bative method for solving QCD and constitutes the ideal framework to describe the structure 
of hadrons in terms of their quark and gluon degrees of freedom. The simple structure of 
the light-front vacuum allows an unambiguous definition of the partonic content of a hadron 
in QCD and of hadronic light-front wavefunctions (LFWFs), the underlying link between 
large distance hadronic states and the constituent degrees of freedom at short distances. 
Thus, one can also solve QCD by diagonalizing the light-front QCD Hamiltonian Hlf- The 
spectrum and light-front Fock-state wavefunctions are obtained from the eigenvalues and 
eigensolutions of the Heisenberg problem Hlf\%1)) = M'^\ip), which becomes an infinite set 
of coupled integral equations for the light-front components ipn = {n\ip) in the Fock expan- 
sion [4, 5]. This nonperturbative method has the advantage that it is frame-independent, 
operates in physical Minkowski space-time, and has no fermion-doubling problem [4]. It has 
been applied successfully in lower space-time dimensions. In practice, however, the resulting 
large matrix diagonalization problem in 3 + 1 space-time has proven to be a daunting task, 
so alternative methods and approximations are necessary. 

The AdS/CFT correspondence between gravity or string theory on a higher- dimensional 



anti-de Sitter (AdS) space and conformal field theories (CFT) in physical space-time [6], 
has led to a semiclassical approximation for strongly-coupled quantum field theories which 
provides physical insights into its nonperturbative dynamics. The correspondence is holo- 
graphic in the sense that it determines a duality between theories in different number of 
space-time dimensions. In practice, the duality provides an effective gravity description in 
a {d+ l)-dimensional AdS space-time in terms of a fiat (i-dimensional conformally-invariant 
quantum field theory defined at the AdS asymptotic boundary [7, 8]. Thus, in principle, one 
can compute physical observables in a strongly coupled gauge theory in terms of a classical 
gravity theory. 

Anti-de Sitter AdS^ space is the maximally symmetric space-time with negative curvature 
and a four-dimensional space-time boundary. The most general group of transformations that 
leave the AdS^+i differential line element 

ds^ = —^ {rj^udx^dx" — dz"^) , (1) 

invariant, the isometry group, has dimensions {d + l){d + 2)/2 (i? is the AdS radius). Five- 
dimensional anti-de Sitter space AdSs has 15 isometries, in agreement with the number of 
generators of the conformal group in four dimensions. Since the AdS metric (1) is invariant 
under a dilatation of all coordinates x^ — )■ Xx^ and z — )■ Xz, it follows that the additional 
dimension, the holographic variable z, acts as a scaling variable in Minkowski space: different 
values of z correspond to different energy scales at which the hadron is examined. As a result, 
a short space-like or time-like invariant interval near the light-cone, x^x'^ — >■ maps to the 
conformal AdS boundary near 2 — )■ 0. This also corresponds to the Q -^ oo ultraviolet (UV) 
zero separation distance. On the other hand, a large invariant four- dimensional interval of 
confinement dimensions x^x^ ~ I/Aq^j-, maps to the large infrared (IR) region of AdS space 
z ~ 1/Aqcd- 

QCD is fundamentally different from conformal theories since its scale invariance is 
broken by quantum effects. A gravity dual to QCD is not known, but the mechanisms 
of confinement can be incorporated in the gauge/gravity correspondence by modifying the 
AdS geometry in the large IR domain z ~ 1/Aqcd, which also sets the scale of the strong 
interactions [9]. In this simplified approach, we consider the propagation of hadronic modes 
in a fixed effective gravitational background asymptotic to AdS space, which encodes salient 
properties of the QCD dual theory, such as the UV conformal limit at the AdS boundary, 
as well as modifications of the background geometry in the large-^ IR region to describe 
confinement. 

The physical states in AdS space are represented by normalizable modes ^p{x,z) = 



e~^^'^^{z), with plane waves along Minkowski coordinates x^ and a profile function $(z) 
along the holographic coordinate z. The hadronic invariant mass states P^P^ = A4.^ are 
found by solving the eigenvalue problem for the AdS wave equation. The modified theory 
generates the point-like hard behavior expected from QCD, instead of the soft behavior 
characteristic of extended objects [9]. It is rather remarkable that the QCD dimensional 
counting rules [10, 11] are also a key feature of nonperturbative models [9] based on the 
gauge/gravity duality. Although the mechanisms are different, both the perturbative QCD 
and the AdS/QCD approaches depend on the leading-twist (dimension minus spin) interpo- 
lating operators of the hadrons and their structure at short distances. 

The gauge/gravity duality leads to a simple analytical and phenomenologically com- 
pelling nonperturbative frame-independent first approximation to the light-front Hamil- 
tonian problem for QCD - "Light-Front Holography" [12]. Incorporating the AdS/CFT 
correspondence [6] as a useful guide, light-front holographic methods were originally intro- 
duced [13, 14] by mapping the Polchinski-Strassler formula for the electromagnetic (EM) 
form factors in AdS space [15] to the corresponding Drell-Yan-West expression at fixed 
light-front time in physical space-time [16, 17]. It was also shown that one obtains identical 
light-front holographic mapping for the gravitational form factor [18] - the matrix elements 
of the energy-momentum tensor, by perturbing the AdS metric (1) around its static solu- 
tion [19]. In the usual "bottom-up" approach to the gauge/gravity duality [20, 21], fields 
in the bulk geometry are introduced to match the chiral symmetries of QCD and axial and 
vector currents become the primary entities as in effective chiral theory. In contrast, in 
light-front holography a direct connection with the internal constituent structure of hadrons 
is established using light-front quantization [12, 13, 14, 18, 22]. 

The identification of higher dimensional AdS space with partonic physics in physical 
space-time is specific to the light front: the Polchinski-Strassler formula for computing tran- 
sition matrix elements is a simple overlap of AdS amplitudes, which maps to a convolution 
of frame-independent light-front wavefunctions. This AdS convolution formula cannot be 
mapped to current matrix elements at ordinary fixed time t, since the instant-time wave- 
functions must be boosted away from the hadron's rest frame - an intractable dynamical 
problem. In fact, the boost of a composite system at fixed time t is only known at weak 
binding. Moreover, the form factors in instant time also require computing the contributions 
of currents which arise from the vacuum in the initial state and which connect to the hadron 
in the final state. Thus instant form wavefunctions alone are not sufficient to compute co- 
variant current matrix elements in the instant form. There is no analog of such contributions 
in AdS. In contrast, there are no vacuum contributions in the light-front formulae for current 



matrix elements - in agreement with the AdS formulae. 

Unlike ordinary instant-time quantization, the Hamiltonian equation of motion in the 
light-front is frame independent and has a structure similar to eigenmode equations in AdS 
space. This makes a direct connection of QCD with AdS/CFT methods possible. In fact, one 
can also study the AdS/CFT duality and its modifications starting from the LP Hamiltonian 
equation of motion for a relativistic bound-state system Hlf\'^) = M'^l'ip) in physical space- 
time [12], where the QCD light-front Hamiltonian Hlf = Pf,P^' = P+P"-Pi, P^ = P^±P\ 
is constructed from the QCD Lagrangian using the standard methods of quantum field 
theory [4]. To a first semiclassical approximation, where quantum loops and quark masses 
are not included, LF holography leads to a LP Hamiltonian equation which describes the 
bound-state dynamics of light hadrons in terms of an invariant impact kinematical variable 
C which measures the separation of the partons within the hadron at equal light-front time 
T = x"^ = x^ + x^. The transverse coordinate ( is closely related to the invariant mass 
squared of the constituents in the LFWF and its off-shellness in the LF kinetic energy, and 
it is thus the natural variable to characterize the hadronic wavefunction. In fact ( is the 
only variable to appear in the relativistic light-front Schrodinger equations predicted from 
holographic QCD in the limit of zero quark masses. The coordinate z in AdS space is thus 
uniquely identified with a Lorentz-invariant coordinate ( which measures the separation of 
the constituents within a hadron at equal light-front time. The AdS/CFT correspondence 
shows that the holographic coordinate z in AdS space is related inversely to the internal 
relative momentum. In fact, light-front holography makes this identification precise. 

Remarkably, the unmodified AdS equations correspond to the kinetic energy terms of 
the partons inside a hadron, whereas the interaction terms in the QCD Lagrangian build 
confinement and correspond to the truncation of AdS space in an effective dual gravity 
approximation [12]. Thus, all the complexities of the strong interaction dynamics are hidden 
in an effective potential U{(), and the central question - how to derive the effective color- 
confining potential U{() directly from QCD, remains open. To circumvent this obstacle, 
the effective confinement potential can be introduced either with a sharp cut-off in the 
infrared region of AdS space, as in the "hard- wall" model [9], or, more successfully, using 
a "dilaton" background in the holographic coordinate to produce a smooth cutoff at large 
distances as in the "soft-wall" model [23]. Furthermore, one can impose from the onset a 
correct phenomenological confining structure to determine the effective IR warping of AdS 
space, for example, by adjusting the dilaton background to reproduce the observed linear 
Regge behavior of the hadronic mass spectrum A^^ as a function of the excitation quantum 
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numbers [23, 24] '^. By using light-front holographic mapping techniques, one also obtains 
a connection between the mass parameter fiR of the AdS theory with the orbital angular 
momentum of the constituents in the light-front bound-state Hamiltonian equation [12]. 
The identification of orbital angular momentum of the constituents is a key element in our 
description of the internal structure of hadrons using holographic principles, since hadrons 
with the same quark content, but different orbital angular momenta, have different masses. 

In our approach, the holographic mapping is carried out in the strongly coupled regime 
where QCD is almost conformal, corresponding to an infrared fixed-point. A QCD infrared 
fixed point arises since the propagators of the confined quarks and gluons in the loop integrals 
contributing to the /3- function have a maximal wavelength [14, 26]; thus, an infrared fixed 
point appears as a natural consequence of confinement. The decoupling of quantum loops 
in the infrared is analogous to QED dynamics where vacuum polarization corrections to the 
photon propagator decouple ai Q^ ^ 0. Since there is a window where the QCD coupling is 
large and approximately constant, QCD resembles a conformal theory for massless quarks. 
One then uses the isometries of AdSs to represent scale transformations within the conformal 
window. We thus begin with a conformal approximation to QCD to model an effective dual 
gravity description in AdS space. The large-distance non-conformal effects are taken into 
account with the introduction of an effective confinement potential as described above. 

Early attempts to derive effective one-body equations in light-front QCD are described 
in reference [27]. We should also mention previous work by 't Hooft, who obtained the 
spectrum of two-dimensional QCD in the large Nc limit in terms of a Schrodinger equation 
as a function of the parton x- variable [28]. In the scale-invariant limit, this equation is 
equivalent to the equation of motion for a scalar field in AdSs space [29]. In this case, there 
is a mapping between the variable x and the radial coordinate in AdSs. 

2 Light-front bound-state Hamiltonian equation of mo- 
tion and light-front holography 

A key step in the analysis of an atomic system, such as positronium, is the introduction of 
the spherical coordinates r,6,(j) which separates the dynamics of Coulomb binding from the 
kinematical effects of the quantized orbital angular momentum L. The essential dynamics 
of the atom is specified by the radial Schrodinger equation whose eigensolutions ipn,L{f) 



■^Using a mean-field mechanism, an effective harmonic confinement interaction was obtained in Ref. [25] 
in a constituent quark model. 



determine the bound-state wavefunction and eigenspectrum. In our recent work, we have 
shown that there is an analogous invariant hght-front coordinate ( which allows one to 
separate the essential dynamics of quark and gluon binding from the kinematical physics 
of constituent spin and internal orbital angular momentum. The result is a single-variable 
light-front Schrodinger equation for QCD which determines the eigenspectrum and the light- 
front wavefunctions of hadrons for general spin and orbital angular momentum [12], thus 
providing a description of the internal dynamics of hadronic states in terms of their massless 
constituents at the same LF time t = x~^ = x^ + x^, the time marked by the front of a light 
wave [3], instead of the ordinary instant time t = x^. 

2.1 Light-front quantization of QCD 

Our starting point is the SU{3)c invariant Lagrangian of QCD 

/:qcd = ^ {t^D, -m)ij- iG^,G"^^ (2) 

where D^ = d^- igsA'^^T" and G"^^ = d^A^ - d^A^ + ^,c'^'"=A|^A;;, with [T'^^T^] = ic^'^^T" 
and a, b, c are SU{3)c color indices. 

One can express the hadron four-momentum generator P = (P"^, P~, P±), P^ = P^±P^, 
in terms of the dynamical fields, the Dirac field ip+, where il!± = A±ip, A± = 7*^7^, and the 
transverse field A_l in the A'^ = gauge [4] quantized on the light-front at fixed light-cone 
time x~^, x^ = x^ ± x^ 

P~ = - dx~cPx±ilj:'j~^ — ?/^_|_ + (interactions), (3) 

2 J io^ 

P+ = /"dx"dW^+7+^9+^+, (4) 

where the integrals are over the null plane r = x"*" = 0, the hyper-plane tangent to the light 
cone. This is the initial-value surface for the fields where the commutation relations are 
fixed. The LF Hamiltonian P~ generates LF time translations 

[V^+(a;),p-] =z^V^+(a;), (6) 

to evolve the initial conditions to all space-time, whereas the LF longitudinal P^ and trans- 
verse momentum P^ are kinematical generators. For simplicity we have omitted from (3-5) 
the contributions from the gluon field Aj^. 



According to Dirac's classification of tlie forms of relativistic dynamics [3], the funda- 
mental generators of the Poincare group can be separated into kinematical and dynamical 
generators. The kinematical generators act along the initial surface and leave the light-front 
plane invariant: they are thus independent of the dynamics and therefore contain no inter- 
actions. The dynamical generators change the light-front position and depend consequently 
on the interactions. In addition to P^ and P^, the kinematical generators in the light-front 
frame are the z-component of the angular momentum J^ and the boost operator K. In 
addition to the Hamiltonian P~, J^ and J^ are also dynamical generators. The light-front 
frame has the maximal number of kinematical generators. 

2.2 A semiclassical approximation to QCD 

A physical hadron in four-dimensional Minkowski space has four-momentum P^ and 
invariant hadronic mass states P^P^ = M.^ determined by the Lorentz-invariant Hamiltonian 
equation for the relativistic bound-state system 

H^p\i,{P)) = M'\i,{P)), (7) 

with Hlp = P^P^ = P^P^ — ?"]_, where the hadronic state \tp) is an expansion in multi- 
particle Fock eigenstates \n) of the free light-front Hamiltonian: {tp) = Xln^^l^)- '^^^ state 
\ip{P^, P_L, J^)) is an eigenstate of the total momentum P+ and P_l and the total spin J^. 
Quark and gluons appear from the light-front quantization of the excitations of the dynam- 
ical fields ip^ and Aj^, expanded in terms of creation and annihilation operators at fixed LF 
time r. The Fock components '?/;„(xj, k_Lj, Aj) are independent of P~^ and P_l and depend 
only on relative partonic coordinates: the momentum fraction Xi = kf / P^ , the transverse 
momentum k^j and spin component Af. Momentum conservation requires Yll=i^i — ^ ^"^^ 
"^^^=1 k±i = 0. The LFWFs -j/^^ provide a frame-independent representation of a hadron 
which relates its quark and gluon degrees of freedom to their asymptotic hadronic state. 
Since for each constituent kf = ykf+mf + k^ > there are no contributions from the 
vacuum. Thus, apart from possible zero modes, the light-front QCD vacuum is the trivial 
vacuum. The constituent spin and orbital angular momentum properties of the hadrons are 
also encoded in the LFWFs. Actually, the definition of quark and gluon angular momen- 
tum is unambiguous in Dirac's front form in light-cone gauge A'^ = 0, and the gluons have 
physical polarization Sg = ±1. 

One can also derive light-front holography using a first semiclassical approximation to 
transform the fixed light-front time bound-state Hamiltonian equation of motion in QCD 
(7) to a corresponding wave equation in AdS space [12]. To this end we expand the initial 
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and final hadronic states in terms of its Fock components. Tlie computation is simplified in 
the frame P = {P+,M^/P+,0±) wliere P^ = P+P~. We find 

// k^ + TTi^ \ 
[dx,] [d^\i^.] Yl [^ -) \Mxi^^±^)\' + (interactions), (8) 

n q '^ 

plus similar terms for antiquaries and gluons {rrig = 0). The integrals in (8) are over the 
internal coordinates of the n constituents for each Fock state 

with phase space norinahzation 

Y^ I [dx^] [d'']^L^] |V^n(x„k^.)|' = 1. (10) 

Each constituent of the light-front wavefunction 'i/'„(xj, k_Lj, Aj) of a hadron is on its 
respective mass shell kf = kfk~ — k^^ = mf, i = 1,2 ■■ -n. Thus k~ = (k^^ + mf) /xiP^ . 
However, the light-front wavefunction represents a state which is off the light-front energy 
shell: P~ — ^" k~ < 0, for a stable hadron. Scaling out P^ = ^^ kf, the invariant mass of 
the constituents A^„ is 

i=l i * 

The functional dependence for a given Fock state is expressed in terms of the invariant mass, 
the measure of the off-energy shell of the bound state of the n-parton LFWF: Ai^— M."^. 

The LFWF ipni^ii k_Lj, Aj) can be expanded in terms of n — 1 independent position coor- 
dinates h±j, j = 1,2, . . . ,n — l, conjugate to the relative coordinates k±i, with ^"^^ b_Lj = 0. 
We can also express Eq. (8) in terms of the internal impact coordinates h±j with the result 



n-l 



Ai'^ = ^^11 / dxj d^h±j 'ijj^{xj,h±j) y^ I — 1 'ipnixj,h±j) + (interactions). 

n J = l>^ q \ ^1 J 

(12) 
The normalization is defined by 



n-l 



lb i n 

YW j dxjd%±j \i^n{xj, h±j)f = 1. (13) 



n j=l 



If we want to simplify further the description of the multiple parton system and reduce 
its dynamics to a single variable problem, we must take the limit of quark masses to zero. 
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Indeed, the underlying classical QCD Lagrangian with massless quarks is scale and conformal 
invariant [30], and consequently only in this limit it is possible to map the equations of motion 
and transition matrix elements to their correspondent conformal AdS expressions. 

To simplify the discussion we will consider a two-parton hadronic bound state. In the 
limit of zero quark mass mg — )■ 



/■I (^x f 

M"^ = / —. r d%±ij*(x,hj_) (-Vl ) V'(a;,b^) + (interactions). 

Jq x{1 -x) J 



(14) 



9 k^ 

For n = 2, the invariant mass is M.n=2 — xii-x) - Similarly, in impact space the relevant 
variable for a two-parton state is C,"^ = x{l — x)b5^. Thus, to first approximation LF dynamics 
depend only on the boost invariant variable A^„ or (, and hadronic properties are encoded 
in the hadronic mode (f){Q from the relation 

^(x,C,¥^) = e^^^X(x)^, (15) 

thus factoring out the angular dependence ip and the longitudinal, X{x), and transverse 
mode </'(C)- This is a natural factorization in the light front since the corresponding canonical 
generators, the longitudinal and transverse generators P^ and P^ and the ^-component of 
the orbital angular momentum J^ are kinematical generators which commute with the LF 
Hamiltonian generator P^ . We choose the normalization (0|0) = j d( |(C|0)P = Pqg^ where 
Pqq is the probability of finding the qq component in the pion light-front wavefunction. The 
longitudinal mode is thus normalized as L ,-, \ = 1. 

<=' JO x{l—x) 

We can write the Laplacian operator in (14) in circular cylindrical coordinates {(, ip) 



I d f ^d\ Id 



^2 



and factor out the angular dependence of the modes in terms of the 5*0(2) Casimir repre- 
sentation L^ of orbital angular momentum in the transverse plane. Using (15) we find [12] 

M' = Jdc<p*{oVc (-^ " ^^ + f ) ^ + ldc<p*{cMom, (17) 

where L = \L^\. In writing the above equation we have summed the complexity of the 
interaction terms in the QCD Lagrangian by the introduction of the effective potential ^(C), 
which is modeled to enforce confinement at some IR scale. The LF eigenvalue equation 
Pfj^P^\(j)) = A^^|0) is thus a light-front wave equation for 

^2 1 _ 4 r 2 \ 

-^ - -^^ + UiO) 0(C) = -MV(C), (18) 
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a relativistic single-variable LF Schrodinger equation. Its eigenmodes (j){() determine the 
hadronic mass spectrum and represent the probability amplitude to find n-partons at trans- 
verse impact separation (, the invariant separation between pointlike constituents within the 
hadron [13] at equal LF time. Thus the effective interaction potential is instantaneous in LF 
time r, not instantaneous in ordinary time t. The LF potential thus satisfies causality, unlike 
the instantaneous Coulomb interaction. Extension of the results to arbitrary n follows from 
the x-weighted definition of the transverse impact variable of the n—1 spectator system [13] 



n-l 






where x = Xn is the longitudinal momentum fraction of the active quark. One can also 
generalize the equations to allow for the kinetic energy of massive quarks using Eqs. (8) 
or (12) [31]. In this case, however, the longitudinal mode X{x) does not decouple from the 
effective LF bound-state equations. 

2.3 Higher spin hadronic modes in AdS space 

We now turn to the formulation of bound-state equations for mesons of arbitrary spin J 
in AdS space '^. As we shall show in the next section, there is a remarkable correspondence 
between the equations of motion in AdS space and the Hamiltonian equation for the rela- 
tivistic bound-state system for the corresponding angular momentum in light-front theory. 

The description of higher spin modes in AdS space is a notoriously difficult problem [34, 
35, 36]. A spin-J field in AdS^+i is represented by a rank J tensor field ^{x^)mi-Mj, which 
is totally symmetric in all its indices. Such a tensor contains also lower spins, which can be 
eliminated by imposing gauge conditions. The action for a spin-J field in AdS^+i space-time 
in presence of a dilaton background field f{z) (the string frame) is given by 



where M, A^ = 1, ■ ■ ■ , d+1, yfg = {R/ zY^^ and Dm is the covariant derivative which includes 
parallel transport. The omitted terms in (20) refer to terms with different contractions. The 



''This section is based on our collaboration with Hans Guenter Dosch. A detailed discussion of higher 
integer and half- integer spin wave equations in modified AdS spaces will be given in Ref. [32]. See also the 
discussion in Ref. [33]. 
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coordinates of AdS are the Minkowski coordinates x'^ and the holographic variable z labeled 
x^'^ = {x^, z). The d + 1 dimensional mass /i is not a physical observable and is a priory an 
arbitrary parameter. The dilaton background field Lp(z) in (20) introduces an energy scale 
in the five-dimensional AdS action, thus breaking its conformal invariance. It is a function 
of the holographic coordinate z which vanishes in the conformal ultraviolet limit z ^ 0. In 
the hard wall model (f = and the conformality is broken by the IR boundary conditions at 

Z = ZQr^ 1/AqCD- 

A physical hadron has plane-wave solutions and polarization indices M along the 3 + 1 
physical coordinates 

^Pix,z),,...,, = e'^^-^^z),,...,,, (21) 

with four- momentum P^ and invariant hadronic mass P^P^'^ = M.^. All other components 
vanish identically: ^z^l2■■■^l,J = ^fi.iz--^j.j = ■ ■ ■ = ^^j^^j...^ = 0. One can then construct an 
effective action in terms of high spin modes $j = $^i^2-/ij5 with only physical degrees of 
freedom [23]. In this case the system of coupled differential equations which follow from (20) 
reduce to a homogeneous equation in terms of the physical field $j. 
In terms of fields with tangent indices 



Mi^A-h Mj ^ _ ( Z 



J 



^AxA2-Aj - (^Ai^^aI ■ ■ -(^Aj ^MiM2-Mj - i^^j ^AiAi-Aj, (22) 

we find the effective action [32] ($j = $^j...^^) 

S = i Jd'^x dz v/^e'^(^) l^g'^'^'dN^jdN'^j - /i^l.^ j , (23) 

containing only the physical degrees of freedom and usual derivatives. Thus, the effect of the 
covariant derivatives in the effective action for spin-J fields with polarization components 
along the physical coordinates is a shift in the AdS mass /i. The vielbein e^^ is defined by 
Qain = €."^^^Vab, where A,B = I,--- ,d + 1 are tangent AdS space indices and rjAB has 
diagonal components (1,— I,--- ,—1). In AdS the vielbein is e^j = {R/z)S^. 

In terms of the AdS field $j = •I'^j...^^ we can express the effective action (23) 

S = ^Jd''xdz^e'P^'^(^g'''''dN^jdN'^j-li'^i), (24) 

where we have defined an effective metric determinant 

^ = iR/zf^'-'', (25) 
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and rescaled the AdS mass fi in (23). Variation of the higher-dimensional action (24) gives 
the AdS wave equation for the spin- J mode $j 

^d-l-2J / -en/ _\ \ / . n\ 2 



oip{z) 



-a 



yd-l-2J 



d. 



fiR\- 



$U 



M^<l>(z) 



where the eigenmode $j is normahzed according to 

dz 



R 



d-l-2J 



1^1^ ^''-'^'A^) 



1. 



The AdS mass is fj, obeys the relation 

{flRf = {T-J){T-d+J), 



(26) 



(27) 



(2^ 



which follows from the scaling behavior of the tangent AdS field near 2 — t- 0, $j ~ 2;^. 

We can also derive (26) by shifting dimensions for a J-spin mode [12, 37]. To this end, we 
start with the scalar wave equation which follows from the variation of (20) for J = 0. This 
case is particularly simple as the covariant derivative of a scalar field is the usual derivative. 
We obtain the eigenvalue equation 

d—l / -infz) \ / n\ 2 



^ifi{z) 



a 



M^) 



yd-l 



a 



/ii?Y 



<^ = M^<^. 



(29) 



A physical spin-J mode $^^...^j with all indices along 3+1 is constructed by shifting dimen- 
sions ^j{z) = {z/R)~-^^{z). It is simple to show that the shifted field $^i/i2---Mj obeys the 
wave equation (26) which follows from (29) upon mass rescaling {^RY -^ if^R)"^ ~ J{d — 
J) + Jzip'{z). 



2.3.1 Non-conformal warped metrics 

In the Einstein frame conformal invariance is broken by the introduction of an additional 
warp factor in the AdS metric in order to include confinement forces 



ds' 



The action is 



S=\jd'xdz^^{gr'9E 



NN' MiM[ 



EJMN 



dx^'dx"" 



{9 
d2 

— e^(") [ri^^dx^'dx" - dz^) . 



(30) 



■9e 



'^'■''^'DN^M....M,DN'^f 



' N^ Kh---MjJ^ N'^ M[---M'j 



2 M^M[ 

/^ 9e 



■9e 



MjM'j 



$ 



$, 



Mi-Mj^M[-M'j 



+ 



, (31) 
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where g^^ = (gE)^^^ and {gE)MN = ^s^^^^Vmn- The flat metric r]MN has diagonal com- 
ponents (1, —1, ■ • • , —1). 

The use of warped metrics is useful to visualize the overall confinement behavior as we 
follow an object in warped AdS space as it falls to the infrared region by the effects of gravity. 
The gravitational potential energy for an object of mass m in general relativity is given in 
terms of the time-time component of the metric tensor (700 

X{z)/2 



e 

z 



V = mcWi9E)oo = mc^R , (32) 



thus we may expect a potential that has a minimum at the hadronic scale zq and grows fast 
for larger values of z to confine effectively a particle in a hadron within distances z ^ zq. In 
fact, according to Sonnenscheim [38] a background dual to a confining theory should satisfy 
the conditions for the metric component g^Q 

dzigoo)\z=zo = 0, 5-00^=^0 7^ 0, (33) 

to display the Wilson loop area law for confinement of strings. 

To relate the results in the Einstein frame where hadronic modes propagate in the non- 
conformal warped metrics (30) to the results in the String- Jordan frame (20), we scale 
away the dilaton profile by a redefinition of the fields in the action. This corresponds to 
the multiplication of the metric determinant ^/gE = (7 ) e*^"'"'"^-*^*^^-'/^ by the contravariant 
tensor (^e)^^^. Thus the result [32] ip{z) -^ ^\{z), or (^ ^ |A for AdSg. 

2.3.2 Effective confining potentials in AdS 

For some applications it is convenient to scale away the dilaton factor in the action by 
a field redefinition [39]. For example, for a scalar field we can shift $ —^ e"'^^'^^, and the 
bilinear component in the action is transformed into the equivalent problem of a free kinetic 
part plus an effective confining potential V{z) which breaks the conformal invariance. ^ For 
the spin-J effective action (24) we find upon the field redefinition $j — )■ e~'^/^<l>j 

S = l jd^xdz v/^(^^^'9^$ja^'$j - ^l^^] - V{z)^^^ 



^In fact, for fcrmions the conformality cannot be broken by the introduction of a dilaton background or 
by exphcitly deforming the AdS metric as discussed above, since the additional warp factor is scaled away by 
a field redefinition. In this case the breaking of the conformal invariance and the generation of the fermion 
spectrum can only be accomplished by the introduction of an effective potential. This is further discussed 
in Sec. 5.2. 
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with effective metric determinant (25) ^Jg^ = {R/zY'^^ ^■^ and effective potential V{z) = 
^U{z), where 

Uiz) = \^"{z) + \^\zf + ^:^^^'(.). (35) 

The action (24) is thus equivalent, modulo a surface term, to the action (34) written 
in terms of the rotated fields $j — )■ e~'''/^$j. The result (35) is identical to the result 
obtained in Ref. [33]. As we will show in the following section, the effective potential (35), 
for d = 4, is precisely the effective light-front potential which appears in Eq. (18), where 
the LF transverse impact variable C, is identified with the holographic variable z. A different 
approach is discussed in Ref. [40] where the infrared physics is introduced by a back-reaction 
model to the AdS metric. See also Refs. [41, 42, 43, 44]. 



2.4 Light-front holographic mapping 

The structure of the QCD light-front Hamiltonian equation (7) for the state \il){P)) is 
similar to the structure of the wave equation (26) for the J- mode $^^...^j in AdS space; 
they are both frame- independent and have identical eigenvalues A^^, the mass spectrum of 
the color-singlet states of QCD. This provides the basis for a profound connection between 
physical QCD formulated in the light-front and the physics of hadronic modes in AdS space. 
However, important differences are also apparent: Eq. (7) is a linear quantum-mechanical 
equation of states in Hilbert space, whereas Eq. (26) is a classical gravity equation; its 
solutions describe spin-J modes propagating in a higher dimensional warped space. Phys- 
ical hadrons are composite, and thus inexorably endowed of orbital angular momentum. 
Thus, the identification of orbital angular momentum is of primary interest in establish- 
ing a connection between both approaches. In fact, to a first semiclassical approximation, 
light-front QCD is formally equivalent to the equations of motion on a fixed gravitational 
background [12] asymptotic to AdSs, where the prominent properties of confinement are 
encoded in a dilaton profile ^{z). 

As shown in Sect. 2.2, one can indeed systematically reduce the LF Hamiltonian eigen- 
value Eq. (7) to an effective relativistic wave equation (18), analogous to the AdS equations, 
by observing that each n-particle Fock state has an essential dependence on the invariant 
mass of the system and thus, to a first approximation, LF dynamics depend only on A^^. 
In impact space the relevant variable is the boost invariant variable C, which measures the 
separation of quarks and gluons, and which also allows one to separate the bound state 
dynamics of the constituents from the kinematics of their internal angular momentum. 
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Upon the substitution z— >C and 

0j(C) = (C/i?)"'/'"-^e^(^)/2<l>j(C), (36) 

in (26), we find for d = A tlie QCD liglit-front wave equation (18) witli effective potential [37] 

uiO = lv"{0 + ^^'(C)^ + ^^'(0, (3T) 

provided that the fifth dimensional mass /i is related to the internal orbital angular momen- 
tum L = max\L^\ and the total angular momentum J^ = L^ + S^ of the hadron. Light-front 
holographic mapping thus implies that the fifth dimensional AdS mass // is not a free pa- 
rameter but scales as 

{fiRf = -(2 - Jf + L^. (38) 

The angular momentum projections in the light-front z direction L^,S^ and J^ are kine- 
matical generators in the front form, so they are the natural quantum numbers to label the 
eigenstates of light-front physics. In general, a hadronic eigenstate with spin J^ in the front 
form corresponds to an eigenstate of J^ = j{j + 1) in the rest frame in the conventional 
instant form. It thus has 2j + 1 degenerate states with J^ = —j, —j + l,---j — l,+j [4], 
thus J represents the maximum value of | J^|, J = max\J^\. 

If L^ < 0, the LF Hamiltonian defined in Eq. (7) is unbounded from below {(J)\Hlp\(J)) < 
and the spectrum contains an infinite number of unphysical negative values of Ai^ which 
can be arbitrarily large. As Ai^ increases in absolute value, the particle becomes localized 
within a very small region near C = 0? since the effective potential is conformal at small (. 
For Ai"^ — !> — oo the particle is localized at C = 0, the particle "falls towards the center" [45]. 
The critical value L = corresponds to the lowest possible stable solution, the ground state 
of the light-front Hamiltonian. For J = the five dimensional mass fi is related to the 
orbital momentum of the hadronic bound state by (fiR)'^ = —4 + L^ and thus (fiR)^ > —4. 
The quantum mechanical stability condition L^ > is thus equivalent to the Breitenlohner- 
Preedman stability bound in AdS [46]. The scaling dimensions are 2 + L independent of J, 
in agreement with the twist-scaling dimension of a two-parton bound state in QCD. It is 
important to notice that in the light-front the 5*0(2) Casimir for orbital angular momentum 
L^ is a kinematical quantity, in contrast to the usual 5*0(3) Casimir L{L + 1) from non- 
relativistic physics which is rotational, but not boost invariant. The 5*0(2) Casimir form 
L^ corresponds to the group of rotations in the transverse LF plane. Indeed, the Casimir 
operator for SO{N) ~ S^-^ is L{L + N - 2). 
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3 Mesons in light-front holography 

Considerable progress has recently been achieved in the study of the meson excitation 
spectrum in QCD from discrete lattices which is a first-principles method [47]. In practice, 
lattice gauge theory computations of eigenvalues beyond the ground-state are very challeng- 
ing. Furthermore, states at rest are not classified according to total angular momentum 
J and Jzi but according to the irreducible representation of the lattice, and thus a large 
basis of interpolating operators is required for the extraction of meaningful data [47]. In 
contrast, the semiclassical light-front holographic wave equation (18) obtained in the pre- 
vious section describes relativistic bound states at equal light-front time with a simplicity 
comparable to the Schrodinger equation of atomic physics at equal instant time. It thus 
provides a framework for a first-order analytical exploration of the spectrum of mesons. In 
the limit of zero-quark masses, the light-front wave equation has a geometrical equivalent to 
the equation of motion in a warped AdS space-time. 

3.1 A hard- wall model for mesons 

As the simplest example we consider a truncated model where quarks propagate freely 
in the hadronic interior up to the confinement scale 1/Aqcd- The interaction terms in the 
QCD Lagrangian effectively build confinement, here depicted by a hard wall potential 

if c< 



U{0=\ ^ -,^ — Aqcd^' (39) 

This provides an analog of the MIT bag model [48] where quarks are permanently confined 
inside a finite region of space. In contrast to bag models, boundary conditions are imposed 
on the boost-invariant variable C? not on the bag radius at fixed time. The wave functions 
have support for longitudinal momentum fraction < x < 1. The resulting model is a 
manifestly Lorentz invariant model with confinement at large distances, while incorporating 
conformal behavior at small physical separation. 

The eigenvalues of the LF wave equation (18) for the potential (39) are determined by 
the boundary conditions (j){z = 1/Aqcd) = 0, and are given in terms of the roots of the 
Bessel functions: A^L.fc = /^l.^Aqcd- Light-front eigenmodes 0(C) are normalized according 
to 






dC<p\C) = i, (40) 



and are given by 

0L,fc(C) = 4^^^,VCJL{CPL,k^QCD) . (41) 

'Jl+L\PL,k) 

Individual hadron states can be identified by their interpolating operators, which are 
defined at the 2 — > asymptotic boundary of AdS space, and couple to the AdS field $(x, z) 
(22) at the boundary limit (See Appendix A). The short-distance behavior of a hadronic 
state is characterized by its twist (canonical dimension minus spin) r = A — a, where a is 
the sum over the constituent's spin a = Yl^=i '^i- The twist of a hadron is also equal to the 
number of its constituent partons n. ^ 

Pion interpolating operators are constructed by examining the behavior of bilinear co- 
variants ipTip under charge conjugation and parity transformation. Thus, for example, a 
pion interpolating operator q'y^'j^q creates a state with quantum numbers J'^'" = O""*", and a 
vector meson interpolating operator q'~f^q a state 1 . Likewise the operator ^7/^75^ creates 
a state with 1+"*" quantum numbers, for example the ai(1260) positive parity meson. If we 
include orbital excitations, the pion interpolating operator is C2+L = qi^l^D^i^ ■ ■ ■ Di^yq. 
This is an operator with total internal orbital momentum L = Yl^i ^ii twist r = 2 + L and 
canonical dimension A = 3 + -L. Similarly the vector- meson interpolating operator is given 
by O^^^ = q-i^'D^i^ ■ ■ ■ Di^yq. The scahng of the AdS field $ (22) near z ^ 0, ^{z) ~ z^, is 
precisely the scaling required to match the scaling dimension of the local meson interpolating 
operators. 

We list in Table 1 the confirmed (4-star and 3-star) isospin 1 = 1 mesons states from the 
updated Particle Data Group (PDG) [49], with their assigned internal spin, orbital angular 
momentum and radial quantum numbers. The 1=1 mesons have quark content \ud), 
4||mm — dd) and \du). The 7 = 1 mesons are the tt, b, p and a mesons. We have not listed 
in Table 1 the 7 = mesons which are a mix of uu, dd and ss, thus more complex entities. 
The light 7 = mesons are 77, r]' , h, h' u, (p, f and /'. This list comprises the puzzling 7 = 
scalar /-mesons, which may be interpreted as a superposition of tetra-quark states with a 
qq, L = 1, S = 1, configuration which couple to a J = state [50]. ^ 

The light 7 = 1 orbital meson spectrum is compared in Fig. 1 with the truncated- 
space model for n = 0. The data is from PDG [49]. The predictions for the lower mass 
mesons are in better agreement with data as compared with Ref. [51], where naive conformal 
dimensions were used instead. However the model fails to account for the pion as a chiral 
A^TT = state. The hard-wall model for mesons is degenerate with respect to the orbital 



^To include orbital L-depcndcnce we make the substitution r — > n + L. 

^The interpretation of the 7ri(1400) is not very clear [50] and is not included in Table 1. Likewise we do 
not include the 7ri(1600) in the present analysis. 
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Table 1: Confirmed 1=1 mesons listed by PDG [49]. The labels L, S and n refer to assigned 
internal orbital angular momentum, internal spin and radial quantum number respectively. For a 
qq State P = (-l)^+\ C = (-1)^+^ . 



L S n J^^ Meson State 
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Figure 1: 1=1 light-meson orbital states in the hard wall model for Aqcd 
doscalar mesons (left) and vector mesons (right). 



0.32 GeV: pseu- 



quantum number L, and thus fails to account for the important L = \L^\ = 1 triplet splitting 
shown in Fig. 1 (right); the ao(980), ai(1260) and 02(1320) states, which corresponds to 
J = |J^| = 0,1,2 respectively. Using the asymptotic expansion of the Bessel function for 
large arguments we find that A4 ~ 2n + L, in contrast to the usual Regge dependence 



M^' 



n 



L found experimentally [50]. As a consequence, the radial modes are not well 



described in the truncated-space model. For example the first radial AdS eigenvalue has a 
mass 1.77 GeV, which is too high compared to the mass of the observed first radial excitation 
of the meson, the 7r(1300). The shortcomings of the hard-wall model described in this section 
are evaded in the soft wall model discussed below, where the sharp cutoff is modified. 



3.2 A soft- wall model for mesons 

As we discussed in Sec. 2.4, the conformal metric of AdS space can be modified within 
the gauge/gravity framework to include confinement by the introduction of an additional 
warp factor or, equivalently, with a dilaton background ^{z), which breaks the conformal 
invariance of the theory. A particularly interesting case is a dilaton profile exp (i/t^s;^) of 
either sign, since it leads to linear Regge trajectories [23] and avoids the ambiguities in the 
choice of boundary conditions at the infrared wall. The corresponding modified metric can 
be interpreted in the higher dimensional warped AdS space as a gravitational potential in 
the fifth dimension 



V{z) 



mc 



mc R- 



(42) 
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Figure 2: Light-front wavefunctions (pn.LiC) iii physical space-time corresponding to a dilaton 
profile exp(Ac^z^): (left) orbital modes (n = 0) and (right) radial modes (L = 0). 

In the case of the negative solution, the potential decreases monotonically, and thus an 
object located in the boundary of AdS space will fall to infinitely large values of z. This is 
illustrated in detail by Klebanov and Maldacena in Ref. [52]. For the positive solution, the 
potential is nonmonotonic and has an absolute minimum at zq ~ l/n. Furthermore, for large 
values of z the gravitational potential increases exponentially, thus confining any object to 
distances {z) ~ 1/k [53, 54]. 

From (37) we obtain for the positive sign confining solution ip = exp (k^z^) the effective 
potential [54] 

U{0 = kX' + 2k\J-1), (43) 

which corresponds to a transverse oscillator in the light-front. For the effective potential 
(43) equation (18) has eigenfunctions 



^nAO 



A+L 



2n! 



{n + L) 



■c 



1/2+L, 



-'<-'/' L^^i^r) 



(44) 



and eigenvalues 



mIj,l 



Ak' 



J + L 



n 



(45) 



The meson spectrum (45) has a string-theory Regge form M.^ r^ n + L: the square of 
the eigenmasses are linear in both the angular momentum L and radial quantum number 
n, where n counts the number of nodes of the wavefunction in the radial variable (. The 



* Similar results arc found in Ref. [33]. 
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LFWFs (44) for different orbital and radial excitations are depicted in Fig. 2. Constituent 
quark and antiquark separate from each other as the orbital and radial quantum numbers 
increase. The number of nodes in the light-front wave function depicted in Fig. 2 (right) 
correspond to the radial excitation quantum number n. 




Figure 3: / = 1 parent and daughter Regge trajectories for the vr-meson family (left) with k 
GeV; and the p-meson family (right) with k = 0.54 GeV. 



0.59 



For the J = L + S meson families Eq. (45) becomes 



MlrS 



n 



T S 

^+2 



(46) 



The lowest possible solution for ra = J = has eigenvalue A^^ = 0. This is a chiral symmetric 
bound state of two massless quarks with scaling dimension 2 and size (C^) ~ l/ft;^, which 
we identify with the lowest state, the pion. Thus one can compute the corresponding Regge 
families by simply adding 4/t;^ for a unit change in the radial quantum number, Ahi^ for a 
change in one unit in the orbital quantum number and 2n'^ for a change of one unit of spin to 
the ground state value of A^^. The spectral predictions for the J = L + S light pseudoscalar 
and vector meson states, listed in Table. 1, are compared with experimental data in Fig. 3 
for the positive sign dilaton model discussed here. The data is from PDG [49]. 

It is important to notice that in contrast to the hard-wall model, the soft-wall model 
with positive dilaton accounts for the mass pattern observed in radial excitations, as well as 
for the triplet splitting for the L = l, J = 0,l,2, vector meson a-states. Using the spectral 
formula (45) we find 

-A^a2(1320) > -^ai(1260) > -^ao(980)- (47) 
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The predicted values are 0.76, 1.08 and 1.32 GeV for the masses of the ao(980), ai(1260) 
and 02(1320) vector mesons, compared with the experimental values 0.98, 1.23 and 1.32 GeV 
respectively. The prediction for the mass of the L = 1, n = 1 state ao(1450) is 1.53 GeV, 
compared with the observed value 1.47 GeV. For other calculations of the hadronic spectrum 
in the framework of AdS/QCD, see Refs. [55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 
68, 69, 70, 71, 72]. ^ 

4 Meson form factors 

A form factor in QCD is defined by the transition matrix element of a local quark 
current between hadronic states. The great advantage of the front form - as emphasized by 
Dirac - is that boost operators are kinematic. Unlike the instant form, the boost operators 
in the front form have no interaction terms. The calculation of a current matrix element 
{P + q\J^\P) requires boosting the hadronic eigenstate from \P) to \P + q), a task which 
becomes hopelessly complicated in the instant form. In addition, the virtual photon couples 
to connected currents which arise from the instant form vacuum. 

In AdS space form factors are computed from the overlap integral of normalizable modes 
with boundary currents which propagate in AdS space. The AdS/CFT duality incorporates 
the connection between the twist scaling dimension of the QCD boundary interpolating op- 
erators to the falloff of the normalizable modes in AdS near its conformal boundary. If 
both quantities represent the same physical observable for any value of the transferred mo- 
mentum squared g^, a precise correspondence can be established between the string modes 
$ in AdS space and the light front wavefunctions of hadrons tpn in physical four dimen- 
sional space-time [13]. In fact, Light-Front Holography was originally derived by observing 
the correspondence between matrix elements obtained in AdS/CFT with the corresponding 
formula using the light-front representation [13]. The same results follow from comparing 
the relativistic light-front Hamiltonian equation describing bound states in QCD with the 
wave equations describing the propagation of modes in a warped AdS space, as shown in the 
previous section [12]. 

4.1 Meson electromagnetic form factor 

In the higher dimensional gravity theory, the hadronic transition matrix element corre- 
sponds to the coupling of an external electromagnetic field A^(x, z), for a photon propagat- 



'For recent reviews see, for example, Refs. [73, 74]. 
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ing in AdS space, with the extended field $p(x, z) describing a meson in AdS is [15] 

Jd^x dz^A'^'ix, z)^*p,{x, 2) Vm$p(x, z) ~ (211)^6^ (P'-P- q) e^(P + PTFkM^). 

(48) 
To simpUfy the discussion we will first describe a model with a wall at 2; ~ 1/Aqcd ^ the 
hard wall model - which limits the propagation of the string modes in AdS space beyond 
the IR separation z ~ 1/Aqcd and also sets the gap scale [9]. We recall from Sec. 2.3 that 
the coordinates of AdSs are the Minkowski coordinates x^ and z labeled x*^ = (x^, 2;), with 
M, A^ = 1, ■ ■ ■ 5, and g is the determinant of the metric tensor. The pion has initial and final 
four momentum P and P' respectively and q is the four-momentum transferred to the pion 
by the photon with polarization e^. The expression on the right-hand side of (48) represents 
the space-like QCD electromagnetic transition amplitude in physical space-time 

{P'\J^mP) = {P + PTFM{q')- (49) 

It is the EM matrix element of the quark current J^ = Cgq'-f'^q, and represents a local coupling 
to pointlike constituents. Although the expressions for the transition amplitudes look very 
different, one can show that a precise mapping of the matrix elements can be carried out at 
fixed light-front time [13, 14]. 

The form factor is computed in the light front from the matrix elements of the plus- 
component of the current J~^ in order to avoid coupling to Fock states with different numbers 
of constituents. Expanding the initial and final meson states |?/'m(P'^, P_l)) in terms of Fock 
components, \iPm) = J2n'^n/M\n), we obtain the DYW expression [16, 17] upon the phase 
space integration over the intermediate variables in the g+ = frame: 

-^M(g^) = X] / [(^^i] [d^^-Li] ^ejllj*^/jyj{Xi,k'_^i, Xi)lpn/M{Xi,k±i, Xi), (50) 

n j 

where the phase space factor [dxi\ [d'^k±i] is given by (9) and the variables of the light cone 
Fock components in the final-state are given by k'j_^ = k_|_j + (l— Xj) q_L for a struck constituent 
quark and k'j^^ = k^^ — Xi q_L for each spectator. The formula is exact if the sum is over all 
Fock states n. The form factor can also be conveniently written in impact space as a sum of 
overlap of LFWFs of the j = 1, 2, ■ ■ ■ ,n — 1 spectator constituents [75] 

n— 1 „ 71—1 

^A/(?^) = X^ n / dxjd^^j exp (^?q_L • ^ Xjh^^j |^„/a/(Xj, b_Lj) f , (51) 

n j=l j=l 

corresponding to a change of transverse momentum XjCi± for each of the n — 1 spectators 
with Er=i b±i = 0. 
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For definiteness we shall consider the vr"*" valence Fock state \ud) with charges e„ = | and 
ej = |. For n = 2, there are two terms which contribute to Eq. (51). Exchanging x o 1 — x 
in the second integral we find 



FMq')=2nJ^'^^^JcdCJo\Cq^l^] \^u^/.{xX)\\ (52) 

where ("^ = x{l — x)b^ and F^+(g = 0) = 1. 

We now compare this result with the electromagnetic form factor in AdS space-time. The 
incoming electromagnetic field propagates in AdS according to A^{x'^, z) = e^{q)e~^'^'^V{q'^, z) 
in the gauge A^ = (no physical polarizations along the AdS variable z). The bulk-to- 
boundary propagator V{q^,z) is the solution of the AdS wave equation for A^j_{x^,z) given 
by (Q2 = _q2 > 0) 

V{Q\z) = zQKi{zQ), (53) 

with boundary conditions [15] 

V{Q'' = Q,z) = V{Q\z = {]) = !. (54) 

The propagation of the pion in AdS space is described by a normalizable mode $p(x'^, z) = 
e~^^'^^{z) with invariant mass P^P^ = M.^ and plane waves along Minkowski coordinates 
x^. Extracting the overall factor {2'kY5^ {P' ~ P ~ q) from momentum conservation at the 
vertex which arises from integration over Minkowski variables in (48), we find [15] 



dz 



Fm = R' / ^V{Q\z)^\z), (55) 



where F{Q'^= 0) = 1. Using the integral representation of V{Q'^, z 

rl 
JO 

we write the AdS electromagnetic form-factor as 



ViQ'', z)= I dxJo\ zQ\l ^-^ I , (56) 




F{Q')=R' dx —Jo\zQ\ \^\z). (57) 



To compare with the light-front QCD form factor expression (52) we use the expression 
of the LFWF (15) in the transverse LF plane, where we factor out the longitudinal and trans- 
verse modes (f){C,) and X{x) respectively. If both expressions for the form factor are to be iden- 
tical for arbitrary values of Q, we obtain (f){C) = {(/R)^^^^{Q and X{x) = ^/x{l — x) [13], 
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where we identify the transverse impact LF variable ( with the holographic variable z, 
z ^ C, = \/x{l — x)|b_L|. ^° Thus, in addition of recovering the expression found in Sec. 2.4 
which relates the transverse mode 0(C) in physical space-time to the field $ in AdS space, 
we find a definite expression for the longitudinal LF mode X{x). Identical results follow 
from mapping the matrix elements of the energy-momentum tensor [18]. 



4.2 Elastic form factor with a dressed current 

The results for the elastic form factor described above correspond to a free current 
propagating on AdS space. It is dual to the electromagnetic point-like current in the Drell- 
Yan-West light-front formula [16, 17] for the pion form factor. The DYW formula is an 
exact expression for the form factor. It is written as an infinite sum of an overlap of LF Fock 
components with an arbitrary number of constituents. This allows one to map state-by-state 
to the effective gravity theory in AdS space. However, this mapping has the shortcoming that 
the multiple pole structure of the time-like form factor does not appear in the time-like region 
unless an infinite number of Fock states is included. Furthermore, the moments of the form 
factor at Q^ = diverge term- by-term; for example one obtains an infinite charge radius [76]. 
This could have been expected, as we are dealing with a massless quark approximation. In 
fact, infinite slopes also occur in chiral theories when coupling to a massless pion. 

Alternatively, one can use a truncated basis of states in the LF Fock expansion with a 
limited number of constituents and the nonperturbative pole structure can be generated with 
a dressed EM current as in the Heisenberg picture, i.e., the EM current becomes modified 
as it propagates in an IR deformed AdS space to simulate confinement. The dressed current 
is dual to a hadronic EM current which includes any number of virtual qq components. The 
confined EM current also leads to finite moments at Q^ = 0, as illustrated in Fig. 4 for the 
EM pion form factor. 

We describe briefly below how to compute a form factor for a confined current in AdS 
space using a soft wall example. However, the actual computation of a form factor in AdS 
has several caveats which we will discuss in Sec. 4.4. 

The effective potential corresponding to a dilaton profile exp(±K^2;^) has the form of a 
harmonic oscillator confining potential n^z"^. The normalizable solution for a meson of twist 
r (the number of constituents for a given Fock component) corresponding to the lowest radial 



^"Extension of the results to arbitrary n follows from the x-weighted definition of the transverse impact 
variable of the n — 1 spectator system given by Eq. (19). In general the mapping relates the AdS density 
$^(2;) to an effective LF single particle transverse density [13]. 
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Figure 4: space-like electromagnetic pion form factor F.,^{q^). Continuous line: confined current, 
dashed line free current. Triangles are the data compilation from Baldini [77], boxes are JLAB 
data [78]. 



n = and orbital L = state is given by 
with normalization 



^^ ^--i^^e--^'^'/^ 



Tir-r 



($-[$-)= f^e-'''''^^{zf = Pr 



(58) 



(59) 



where Pr is the probability for the twist r mode (58). This agrees with the fact that the field 
$^ couples to a local hadronic interpolating operator of twist r defined at the asymptotic 
boundary of AdS space (See Appendix A), and thus the scaling dimension of $'^ is r. 
In the case of a soft-wall potential the EM bulk-to-boundary propagator is [14, 79] 



where U{a, b, c) is the Tricomi confluent hypergeometric function 

POO 

T{a)U{a,b,z)= / e-'H^-^l + ty-^-^dt. 



(60) 



(61) 



The modified current VlQ"^, z), Eq. (60), has the same boundary conditions (54) as the free 
current (53), and reduces to (53) in the limit Q^ ^ oo [14]. Eq. (60) can be conveniently 
written in terms of the integral representation [79] 

"^ dx 



ViQ^z) 



2 2 
K Z 






-K^Z^x/(l — x) 



X) 



(62) 
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Substituting in (55) the expression for the hadronic state (58) with twist r and the bulk- 
to-boundary propagator (62), we find that the corresponding elastic form factor for a twist 
r Fock component Fr{Q^) (Q^ = _g2 > q) [^4] 

FAQ') = ^7w ^?x -^ ^r^^ (63) 



M'i I \^ ' M", I V" ' M^T^ 



which is expressed as a r — 1 product of poles along the vector meson Regge radial trajectory. 
For a pion, for example, the lowest Fock state - the valence state - is a twist-2 state, and 
thus the form factor is the well known monopole form [14]. Thus the mean-square charge 
radius of the pion (r^) = Q/Ai^ in the valence approximation. For Aip ^ 770 MeV we find 
(r^r) — 0.63 fm, compared with the experimental value (r^) = 0.672 ± 0.008 fm [49]. In 
contrast, the computation with a free current gives the logarithmically divergent result ^^. 

The remarkable analytical form of (63), expressed in terms of the p vector meson mass and 
its radial excitations, incorporates the correct scaling behavior from the constituent's hard 
scattering with the photon and the mass gap from confinement. 

4.3 Effective wave function from holographic mapping of a con- 
fined current 

It is also possible to find a precise mapping of a confined EM current propagating in a 
warped AdS space to the light-front QCD Drell-Yan-West expression for the form factor. In 
this case we we find an effective LFWF, which corresponds to a superposition of an infinite 
number of Fock states generated by the "dressed" confined current. For the soft- wall model 
this mapping can be done analytically. 

The form factor in light-front QCD can be expressed in terms of an effective single- 
particle density [75] 

F{Q^)= f dxp{x,Q), (65) 

where 

/■oo 

p{x,Q) = 2Tr bdbJoibQ{l-x))\'4;{x,b)\'^, (66) 

Jo 



^"'^The logarithmically divergent result does not appear in the hard- wall model if one uses Neumann bound- 
ary conditions. In this case the EM current is confined and (rj) ~ I/Aq^q. A discussion of the pion form 
factor including chiral symmetry breaking effects in the hard- wall model is given in Refs. [80] and [81]. 
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for a two-parton state (6 = |b_L|). 

We can also compute an effective density on the gravity side corresponding to a twist r 
hadronic mode $,- in a modified AdS space. For the soft-wall model the expression is [14] 



r,2 
\r-2 



p{x,Q) = {T-l){l-xy~'x^. (67) 

To compare (67) with the QCD expression (66) for twist-two we use the integral 

u du Joiau) e-/^"' = ^ e""'/^'^, (68) 

2p 

and the relation x^ = e^^'^^^\ We find the effective two-parton LFWF 



in impact space. The momentum space expression follows from the Fourier transform of (69) 
and it is given by 



an (1) 

K[l — X) 



^^y}^^-^l/2.^il^.)Hni^^)_ (71) 

k(1 — X) 

The effective LFWF encodes nonperturbative dynamical aspects that cannot be learned 
from a term-by-term holographic mapping, unless one includes an infinite number of terms. 
Furthermore, it has the right analytical properties to reproduce the bound state vector meson 
pole in the time-like EM form factor. Unlike the "true" valence LFWF, the effective LFWF, 
which represents a sum of an infinite number of Fock components, is not symmetric in the 
longitudinal variables x and 1 — x for the active and spectator quarks, respectively. 

4.4 Some caveats computing matrix elements in AdS/QCD 

The positive dilaton background exp(-|-K^2^) used in Sec. 3.2 leads to a successful de- 
scription of the meson spectrum in terms of the internal quantum numbers n, L and S, and 
has been preferred for computations in the framework of light-front holography, where the in- 
ternal structure of hadrons is encoded in the wave function. The positive dilaton background 
has been discussed in the literature [33, 53, 54, 82, 83] since it has the expected behavior of a 
model dual to a confining theory [38, 52]. This solution was studied in Ref. [23] but discarded 
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in the same paper, as it leads to a spurious massless scalar mode in the two-point correla- 
tion function for vector mesons [84], and a dilaton field with opposite sign, exp(— k^z^), was 
adopted instead [23]. However, using the results of Sec. 2.3.2, one can readily show that the 
difference in the effective potential U{z) corresponding to positive and negative dilaton fac- 
tors exp{±K,'^z'^) simply amounts to a 2;- independent shift in the light-front effective potential 
U, which in fact vanishes in the vector meson J = 1 channel. From (35) 

AUiz) = U^{z) - U.^{z) = ip"{z) + ^^^:^±i<^'(2), (72) 

At 

in agreement with the results found in Ref. [33]. 

For the dilaton profile ip = k'^z'^ we find for d = A 

AU = A{J-1)k^. (73) 

Therefore, from the point of view of light-front physics, plus and minus dilaton soft-wall 
solutions are equivalent upon a redefinition of the eigenvalues for J 7^ 1. For J = 1 the 
effective potential is U = k^z"^, identical for the plus and minus solutions [39]. Thus, the 
five-dimensional effective AdS action for a conserved EM current Vm in presence of a confining 
potential U = kS^ [39] 

S = Jd'x dz^g Qfm^F*^^ - ^VmV'^ , (74) 

where Fmn = QmVn — OnVm, only differs by a surface term from the action corresponding 
to plus or minus dilaton profiles. Equivalently, one can start from the five- dimensional 
action (74). Upon the field redefinition Vm — ^ e^*^ ^ ^"^Vm one obtains the five-dimensional 
actions corresponding to plus or minus dilaton solutions, which differ from (74) only by a 
surface term. Consequently, essential physics cannot dependent on the particular choice of 
the dilaton sign. 

Another difficulty found in the holographic approach to QCD is that the vector meson 
masses obtained from the spin-1 equation of motion do not match the poles of the dressed 
current when computing a form factor. The discrepancy, in the case of the pion, is an 
overall factor of v2 between the value of the gap scale which follows from the spectrum or 
from the computation of the pion form factor in the valence state approximation.^^ This 



^^This discrepancy is also present in the gap scale if one computes the spectrum and form factors without 
recourse to holographic methods, for example using the semi-classical approximation of Ref. [12]. In this 
case a discrepancy of a factor factor v2 is also found between the spectrum and the computation of space-like 
form factors. 
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is quite puzzling, since the same discrepancy is also found, for example, when computing a 
space-like form factor using the Drell-Yan-West expression, which is an exact expression if 
all Fock states are included. In AdS conserved currents are not renormalized and correspond 
to five dimensional massless fields propagating in AdS according to the relation (/x-R)^ = 
(A— p)(A+p— 4) for ap form. In the usual AdS/QCD framework [20, 21] this corresponds for 
p = ltoA = 3orl, the canonical dimensions of an EM current and the massless gauge field 
respectively. Normally, one uses a hadronic interpolating operator with minimum twist r to 
identify a hadron and to predict the power-law fall-off behavior of its form factors and other 
hard scattering amplitudes [9]; e.g., for a two-parton bound state r = 2. However, in the case 
of a current, one needs to use an effective field operator with dimension A = 3. The apparent 
inconsistency between twist (28) and canonical dimension is removed by noticing that in the 
light-front one chooses to calculate the matrix element of the twist-3 plus component of the 
"good" current J~^ [13, 14], in order to avoid coupling to Fock states with different numbers 
of constituents [16, 17]. 

As described in Sec. 2.4, light front holography provides a precise relation of the fifth- 
dimensional mass /i with the total and orbital angular momentum of a hadron in the trans- 
verse LF plane (/i-R)^ = — (2 — J)^ + L^ (38). Thus the poles computed from the AdS wave 
equations for a conserved current fiR = 0, correspond to a J = L = 1 twist-3 state. Fol- 
lowing this, we can compute the mass of the radial excitations of the twist-3 vector family 
J = L = 1 using Eq. (45). The result is 



Mlj=,^L=i = '^^\n + l), (75) 

which is identical with the results obtained in Ref. [23], since, as explained above, the meson 
spectrum computed with positive or negative dilaton solutions is indistinguishable for J = 1. 
The twist-3 computation of the space-like form factor, involves the current J"*", and the 
poles given by (75) do not correspond to the physical poles of the twist-2 transverse current 
J_L present in the annihilation channel, namely the J = 1,L = state. In this case Eq. (45) 
gives for the twist-2, J = 1, L = vector family the result 

-Mn,j=i,L=o = 4«:' (^ + • ^^6) 

Thus, to compare with physical data one must shift in (63) the twist-2 poles given by (75) 
to their physical positions (76). When the vector meson masses are shifted to their physical 
values the agreement of the predictions with observed data is very good [85]. We presume 
that the problem arises because of the specific truncation used. 
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4.5 Meson transition form factors 

The photon-to-meson transition form factors ^^ (TFFs) Fm^jIQ'^) measured in 77* — )■ M 
reactions have been of intense experimental and theoretical interest. The pion transition 
form factor between a photon and pion measured in the e~e~ — ?■ e~e"7r'' process, with one 
tagged electron, is the simplest bound-state process in QCD. It can be predicted from first 
principles in the asymptotic Q'^ ^ 00 limit [88]. More generally, the pion TFF at large 
Q^ can be calculated at leading twist as a convolution of a perturbative hard scattering 
amplitude Th{'J'j* — )■ qq) and a gauge-invariant meson distribution amplitude (DA), which 
incorporates the nonperturbative dynamics of the QCD bound-state [88]. 

The BaBar Collaboration has reported measurements of the transition form factors from 
7*7 — >■ M process for the vr^ [89], rj, and rj' [90, 91] pseudoscalar mesons for a momentum 
transfer range much larger than previous measurements [92, 93]. Surprisingly, the BaBar 
data for the 7r''-7 TFF exhibit a rapid growth for Q^ > 15 GeV^, which is unexpected from 
QCD predictions. In contrast, the data for the 77-7 and t]'-'-/ TFFs are in agreement with 
previous experiments and closer in agreement with theoretical predictions. Many theoretical 
studies have been devoted to explaining BaBar's experimental results [94, 95, 96, 97, 98, 99, 
100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110]. 

The pion transition form factor Ft^^{Q'^) can be computed from first principles in QCD. 
To leading leading order in asiQ"^) and leading twist the result is [88] (Q^ = —q^ > 0) 



V -t'wjyUj ) /— 



dx- 



{x,xQ) 



X 



l + 0(a„ — 



(77) 



where x is the longitudinal momentum fraction of the quark struck by the virtual photon 
in the hard scattering process and x = 1 — x is the longitudinal momentum fraction of the 
spectator quark. The pion distribution amplitude (/)(x, Q) in the light-front formalism [88] 
is the integral of the valence qq LFWF in light-cone gauge A'^ = 



(x,Q) 



167r3 



i:gq/n{x,kj 



(7J 



and has the asymptotic form [88] (j){x, Q — )■ 00) = -\/3/^x(l — x); thus the leading order QCD 
result for the TFF at the asymptotic limit is obtained [88], 



Q'F^,{Q' ^ 00) = 2U 



(79) 



^■^This section is based on our collaboration with Fu-Guang Cao. Further details are given in [86, 87]. 
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To describe the two-photon processes 77* — )■ M , using hght-front holographic methods 
similar to those described in Sec. 4, we need to explore the mathematical structure of higher- 
dimensional forms in the five dimensional action, since the amplitude (48) can only account 
for the elastic form factor FuiQ"^) [87]. For example, in the five-dimensional AdS action there 
is an additional Chern-Simons (CS) term in addition to the usual Yang-Mills term F'^ [8]. In 
the case of a U{1) gauge theory the CS action is of the form e^^'^^^^ AidM^Ndp^Q- The CS 
action is not gauge invariant: under a gauge transformation it changes by a total derivative 
which gives a surface term. The CS form is the product of three fields at the same point in 
five-dimensional space corresponding to a local interaction. Indeed the five- dimensional CS 
action is responsible for the anomalous coupling of mesons to photons and has been used 
to describe, for example, the u — t- -k^ [HI] decay as well as the 77* — )■ tt" [112, 113] and 
7*p° — 7- 7r° [114] processes. ^^ 

The hadronic matrix element for the anomalous electromagnetic coupling to mesons in 
the higher gravity theory is given by the five-dimensional CS amplitude 

dS I dze'''^''''^ALdMANdpAQ 

~ (27r)^5(4) {P + q-k) F^^(g2)eA«-P-e^(g)p,e^(A;)g,, (80) 

which includes the pion field as well as the external photon fields by identifying the fifth 
component of A with the meson mode in AdS space [116]. In the right-hand side of (80) q 
and k are the momenta of the virtual and on-shell incoming photons respectively with corre- 
sponding polarization vectors e^(q') and e^(/c) for the amplitude 77* — )■ vr^. The momentum 
of the outgoing pion is P. 

We now compare the QCD expression on the right-hand side of (80) with the AdS 
transition amplitude on the left-hand side. As for the elastic form factor discussed in Sec. 
4.1, the incoming off-shell photon is represented by the propagation of the non-normalizable 
electromagnetic solution in AdS space, A^{x^, z) = e^{q)e~^'^'^V{q'^, z), where V^q"^, z) is the 
bulk-to-boundary propagator with boundary conditions (54) V{q^ = 0, z) = V{q^, z = Q) = 
1. Since the incoming photon with momentum k is on its mass shell, A;^ = 0, its wave function 
is A^{x^, z) = e^{k)e^^'^ . Likewise, the propagation of the pion in AdS space is described by 
a normalizable mode ^p{x^,z) = e~*^'^$^(z) with invariant mass P^P^ = M.\ = in the 
chiral limit for massless quarks. The normalizable mode ^Tt{z) scales as $7r(-2) — ^ z^ in the 
limit 2; — )■ 0, since the leading interpolating operator for the pion has twist two. A simple 
dimensional analysis implies that A^ ~ ^^^{z)/ z, matching the twist scaling dimensions: 



^^Thc anomalous EM couplings to mesons in the Sakai and Sugimoto model is described in Ref. [115]. 
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two for the pion and one for the EM field. Substituting in (80) the expression given above 

for the the pion and the EM fields propagating in AdS, and extracting the overall factor 

(27r)^6^ {P' — Q — k) upon integration over Minkowski variables, we find {Q^ = —q^> 0) 

1 /"^ rl7 
F.,m = 7r -^.{^)V{Q\z), (81) 

ZTT Jq z 

where the normalization is fixed by the asymptotic QCD prediction (79). We have defined 
our units such that the AdS radius R= 1. 

Since the LP mapping of (81) to the asymptotic QCD prediction (79) only depends 
on the asymptotic behavior near the boundary of AdS space, the result is independent of 
the particular model used to modify the large z IR region of AdS space. At large enough 
Q, the important contribution to (79) only comes from the region near z ~ 1/Q where 
$(z) = 271 f^z"^ + 0{z'^). Using the integral J^ dxx"'Ki{x) 
recover the asymptotic result (79) 



2°-2a [r(f)] , Re(a) > 1, we 



with the pion decay constant fj, [87] 

fn 



V. 



1 9,$^(z) 



^'^ 



^2) 



(83) 

A simple analytical expression for the pion TFF can be obtained from the "soft-wall" 
holographic model described in Sec. 4.2. Using (58) to describe the twist-two pion valence 
wave function in AdS space we find 



\c^ r q^^jXicj 



4 

7! 



\x^ 







X 



1 — exp 



A-K^flx 



(84) 



where 0(a:) = v3/,r3^(l — 2;) is the asymptotic QCD distribution with /^r the pion decay 
constant and Pqq is the probability for the valence state. Remarkably, the holographic 
result for the pion TFF factor given by (84) for Pqq = 1 is identical to the results for the 
pion TFF obtained with the exponential light-front wave function model of Musatov and 
Radyushkin [117] consistent with the leading order QCD result [88]. Since the pion field is 
identified as the fifth component of Am-, the CS form e^^'^^^'^ AlOmAmOpAq is similar in 
form to an axial current; this correspondence can explain why the resulting pion distribution 
amplitude has the asymptotic form. ^^ 



^^In Ref. [112] the pion TFF was studied in the framework of a CS extended hard- wall AdS/QCD model 
with Az ^ dz^{z). The expression for the TFF which follows from (80) then vanishes at Q^ = 0, and has 
to be corrected by the introduction of a surface term at the IR wall [112]. However, this procedure is only 
possible for a model with a sharp cutoff. 
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Figure 5: The 77* -> 7r° transition form factor F.„^{Q'^) (left) and Q'^F.„^{Q'^) (riglit). The dotted 
curve is the asymptotic result. The dashed and solid curves include the effects of using a confined 
EM current for twist-two and twist-two plus twist-four respectively. The data are from [89, 92, 93]. 



Taking Pqq = 0.5 in (84) one obtains a result in agreement with the Adler, Bell and 
Jackiw anomaly result which agrees within a few percent with the observed value obtained 
from the the decay 7r° -^ 77. This suggests that the contribution from higher Fock states 
vanishes at Q = in this simple holographic confining model. Thus (84) represents a 
description of the pion TFF which encompasses the low-energy nonperturbative and the 
high-energy hard domains, but includes only the asymptotic distribution amplitude of the 
qq component of the pion wave function at all scales. The results from (84) for Pgg = 0.5 
are shown in Fig. 5. Also shown in Fig. 5 are the results for the free current approximation 
(which corresponds to the asymptotic result) with Pgg = 0.5 and a twist-two plus twist-four 
model [87] with Pgg = 0.915, and -P^g^g = 0.085. The calculations [87] agree reasonably well 
with the experimental data at low- and medium-Q^ regions {Q^ < 10 GeV^), but disagree 
with BaBar's large Q^ data. 

The rj and rj' mesons result from the mixing of the neutral states rjs and rji of the SU(3)i7' 
quark model. The TFFs for the 77 and rj' mesons have the same expression as the pion 
transition form factor, except for an overall multiplying factor cp = 1, 4=, and ^^ for the 
7r°, r/g and 771, respectively [87]. The results for the rj and rj' transitions form factors are 
shown in Fig. 6. The calculations agree very well with available experimental data over a 
large range of Q^. The rapid growth of the large Q^ data for the pion-photon transition 
form factor reported by the BaBar Collaboration is difficult to explain within the current 
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plus twist-two plus twist-four respectively. Same for the 77* — )• rj' transition form factor Q'^F^i^{Q'^) 
(right). The data are from [89, 92, 93]. 



framework of QCD. The analysis presented here thus indicates the importance of additional 
measurements of the pion-photon transition form factor at large Q^. 



5 Baryons in light-front holography 

The study of the excitation spectrum of baryons is one of the most challenging aspects 
of particle physics. In fact, dedicated experimental programs are in place to determine the 
spectrum of nucleon excitations and its internal structure. Important computational efforts 
in lattice QCD aim to the reliable extraction of the excited nucleon eigenstates. Lattice 
calculations of the ground state light hadron masses agree with experimental values within 
5% [47]. However, the excitation spectrum of the nucleon represents a formidable challenge 
to lattice QCD due to the enormous computational complexity required beyond the leading 
ground state configuration [118]. Moreover, a large basis of interpolating operators is required 
since excited nucleon states are classified according to irreducible representations of the 
lattice, not the total angular momentum. 

As we shall discuss below, the analytical exploration of the baryon spectrum and nu- 
cleon form factors, using light-front gauge/gravity duality ideas, leads, in contrast, to simple 
formulas and rules which describe quite well the systematics of the established light-baryon 



37 



resonances and elastic and transition nucleon form factors, which can be tested against new 
experimental findings. The gauge/gravity duality can give us important insights into the 
strongly coupled dynamics of nucleons using simple analytical methods. 

We can extended the holographic ideas to spin-| hadrons by considering the propagation 
of spin- 1 Dirac modes in AdS space [15]. The action for a Dirac field in AdS^+i is 

Sf = Jd'xdz^ f^^e'jT^DM^ - '-{DM^)e^T^^ - ^^^A , (85) 

where y/g = (^) and ej( is the inverse vielbein, e^ = (-|) 6^. The covariant derivative 
of the spinor field is D^ = Dm — ^^ff^AB where T^ab are the generators of the Lorentz 
group in the spinor representation, Tjab = \ [^a,'^b], and the tangent space Dirac matrices 
obey the usual anti-commutation relation |r'^,r'^} = r]^^. For d even we can choose the 
set of gamma matrices F^ = (F^,F2) with F^ = — F^ = FqFi ■ ■ ■ F^^^i. For c? = 4 we 
have Ta = (7/^,-^75), where 7^ and 75 are the usual 4-dimensional Dirac matrices with 
75 = ^7o7i7273 and 7I = +1. The spin connection in AdS is wff = {ri^^Sfj — r]^^Sfj) /z, 
thus the equation of motion (ieA^^DM — fi) "^ = leads to the Dirac equation in AdS space 



zv^^'^YMd^ + ^r, ) - f,R 



^ = 0, (86) 



where the d + 1 dimensional mass /i is a priory an arbitrary parameter. ^^ 

One can also take as starting point the construction of light-front wave equations in 
physical space-time for baryons by studying the LF transformation properties of spin 1/2 
states [121]. The light- front wave equation describing baryons is a matrix eigenvalue equation 
DlfH^) = ■Ml'ip) with Hlf = D\p. In a 2 x 2 spinor component representation 

7 J . I 1 

dC C 

-^^_ + ^^_ = Mij+. (87) 

As shown below, we can identify u with the orbital angular momentum L: u = L + 1. 
Upon the substitution z ^ ( and 

^(x, z) = e-'^-''z'^i/j{z)u{P), (88) 



^^The spinor action (85) is often complemented by an additional surface term in the UV boundary [119] 
linic^o / d x.y^'^'^ where g^ is the metric induced in the boundary surface by the metric g of AdS^+i. The 
additional term is required to preserve the 0{d+ 1,1) isometry group of AdS^+i and to compute a two-point 
correlation function in the conformal boundary theory [120]. The equation of motion (86) is not modified 
by the surface term. 
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in (86) we recover for d = 4 its LF expression (87), provided that \fiR\ 



i The baryon 



invariant mass is P^P'^ = Ai^ and the spinor u{P) is a four-dimensional spinor which obeys 



the Dirac equation (^ — A^)m(P) = 0. Thus the eigenvalue equation Hlf4'± 
the upper and lower components leads to the wave equation 



d^ 



l-Au' 



dC 



and 



with solutions 



dC^ 



4(2 
l-4(t/ + l)2 



^+(C) = -M'^+(C), 



M.'^il)± for 



(89) 



V;_(C) = A^'V^-(C), 



^+ ~ VCJ.(C-M), V^^ ~ v^J.+i(C>l). 



(90) 



(91) 



The solution of the spin-| Rarita-Schwinger equation for the field \1/a/ in AdS space 
is more involved, but considerable simplification occurs in the '^ z = ^ gauge for physical 
polarization along Minkowski coordinates \1/^, where it becomes similar to the spin-| solu- 
tion [122, 123]. 



5.1 A hard- wall model for baryons 

The hermiticity of the LF Dirac operator Dlf in the eigenvalue equation Dlf\'4^) = M\ip) 
implies that the surface term ip^{()ip-{() —ip'^^Oip^lO should vanish at the boundary. Thus 
in a truncated space holographic model, the light front modes ip^ or ^__ should vanish at the 
boundary C = a-^d C = Co- This condition fixes the boundary conditions and determine 
the baryon spectrum in the truncated hard-wall model. A similar surface term arises when 
one computes the equation of motion from the action (85). In fact, integrating by parts (85) 
and using the equation of motion we find 



Sf 



— lim 

e-5-O 






^_^_ 



20 



(92) 



where \l/± = 2(1=!= 75) ^) ^^^ R has units R = 1. The baryon mass spectrum thus fol- 
lows from the LF "bag" boundary conditions ip± (Co) = or the AdS boundary conditions 
^± (-2^0) = at the IR value, zq = 1/Aqcd, where the LF invariant impact variable ( (19) is 
identified with the AdS holographic coordinate z, z ^ (. We find 

M+ = I3^,k Aqcd, -M" = l3^+i,k Aqcd, (93) 

with a scale-independent mass ratio determined by the zeros of Bessel functions (3i^k- 
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In the usual AdS/CFT correspondence the baryon is an SU{Nc) singlet bound state 
of Nc quarks in the large Nq limit. Since there are no quarks in this theory, quarks are 
introduced as external sources at the AdS asymptotic boundary [124, 125]. The baryon is 
constructed as an Nc baryon vertex located in the interior of AdS. In this top-down string 
approach baryons are usually described as solitons or Skyrmion-like objects [126, 127]. In 
contrast, the bottom- up light-front holographic approach described here is based on the 
precise mapping of AdS expressions to light-front QCD. Consequently, we construct baryons 
corresponding to Nq = 3 not Nc -^ oo. The corresponding interpolating operator for an 
Nc = 3 physical baryon O^+l = '4'D{i^ ■ ■ ■ Di^tpDi^^^ . . . De^ytp, L = Y,T=i ^i' is a twist-3, 
dimension 9/2 + L with scaling behavior given by its twist-dimension 3 + L. We thus require 
z/ = L + 1 to match the short distance scaling behavior. One can interpret L as the maximal 
value of I L^ I in a given LF Fock state. 

In the case of massless quarks, the nucleon eigenstate {u± = | (1 ± 75) u) 

= C^/C{MCM)u^ + J.+i{CM)u.), (94) 

has components ip+ and ip- with different orbital angular momentum, L^ = and L^ = +1, 
combined with spin components S^ = +1/2 and S^ = —1/2 respectively, but with equal 
probability ^^ 

''dC\MCW= [dC\i^4C)\', (95) 



a manifestation of the chiral invariance of the theory for massless quarks. Thus in light- 
front holography, the spin of the proton is carried by the quark orbital angular momentum: 
J^ = (L^) = ±1/2 since (^ S^) = [128], and not by its gluons. 

An important feature of bound-state relativistic theories is that hadron eigenstates have 
in general Fock components with different L components. In the holographic example dis- 
cussed above, the proton has S and P components with equal probability. In the case of 
QED, the ground state IS* state of the Dirac-Coulomb equation has both L = and L = 1 
components. By convention, in both light-front QCD and QED, one labels the eigenstate 
with its minimum value of L. For example, the symbol L in the light-front AdS/QCD spec- 
tral prediction for mesons (46) refers to the minimum L (which also corresponds to the 
leading twist) and S is the total internal spin of the hadron. 



^'^For the truncated-space model, (95) follows from the identity /^ xdx [j^(x/3) — J^^i{x/3)] ~ 
Ja{P)Ja+i{P)/ P, independently of the component wavefunction chosen to fix the boundary conditions at 

C = Co. 
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Table 2: Classification of confirmed baryons listed by the PDG [49]. The labels L, S and n refer 
to the internal orbital angular momentum, internal spin and radial quantum number respectively. 
The even-parity baryons correspond to the 56 multiplet of SU{6) and the odd-parity to the 70. 



L 


S 


n 


Baryon State 





1 

2 





Arl+(940) 





1 
2 


1 


iVi+(1440) 





1 
2 


2 


Ari+(1710) 





3 

2 





A|+(1232) 





3 

2 


1 


A|+(1600) 


1 


1 
2 





iVi"(1535) iV|"(1520) 


1 


3 

2 





iV|"(1650) iV|^(1700) iV|^(1675) 


1 


1 
2 





Ai"(1620) Af~(1700) 


2 


1 
2 





iV|+(1720) Ar|+(1680) 


2 


1 
2 


1 


Ar|+(1900) 


2 


3 

2 





Ai+(1910) A|+(1920) A|+(1905) A|+(1950) 


3 


1 
2 





Nf Nl' 


3 


3 

2 





iVf" N\~ iV|"(2190) iV|"(2250) 


3 


1 
2 





Af- ^- 


4 


1 
2 





A^|+ iV|+(2220) 


4 


3 

2 





A|+ A|+ Af A^+(2420) 


5 


1 
2 





iVf- N^ 


5 


3 

2 





Nl~ iVf" iV^~(2600) Nf~ 
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We list in Table 2 the confirmed (3-star and 4-star) baryon states from the updated 
Particle Data Group [49]. ^^ To determine the internal spin, internal orbital angular mo- 
mentum and radial quantum number assignment of the A^ and A excitation spectrum from 
the total angular momentum-parity PDG assignment, it is convenient to use the conven- 
tional SU{6) D SU{3) flavor X SU{2)spin multiplet structure, but other model choices are also 
possible [130]. ^^ 
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Figure 7: Light baryon orbital spectrum (n = 0) for Aqcd = 0.25 GeV. Predictions for the nucleons 
(left) and for the A trajectories (right). 



We show in Fig. 7 the model predictions for the orbital excitation spectrum of baryons 
which follows from the boundary conditions il)± (C = 1/Aqcd) = in a truncated-space model 
in the infrared region [51]. '^^ The figure shows the predicted orbital spectrum of the nucleon 
and A orbital resonances for n = 0. The only parameter is the value of Aqqd which we take as 
0.25 GeV. Orbital excitations are approximately aligned along two trajectories corresponding 
to even and odd parity states, with exception of the A| (1620) and A| (1700) states which 
are in the same trajectory. The spectrum shows a clustering of states with the same orbital 
L, consistent with a strongly suppressed spin-orbit force. This remarkable prediction for 



^®A recent exploration of the properties of baryon resonances derived from a multichannel partial wave 
analysis [129] report additional resonances not included in the Review of Particle Properties [49]. 

^^In particular the A| (1930) state (not shown in Table 2) has been given the non-SU{6) assignment 
S = 3/2, L = 1, n = 1 in Ref. [130]. This assignment will be further discussed in the section below. 

^°The results shown here in Fig. 7 give better results for the lower mass baryons as compared with Ref. 
[51] where naive conformal dimensions were used instead. 
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the baryons is not a peculiarity of the hard-wall model, but is an important property of 
light-front holographic models. 

In the quark-diquark model of Jaffe and Wilczek [131], nucleon states with S* = 1/2 in 
Fig. 7 (a) correspond to "good" diquarks, S = 3/2 nucleons and all the A states in Fig. 
7 (b) to "bad" diquarks, with exception of the A(1930) which does not follow the simple 
3g quark-diquark pattern. As for the case for mesons discussed in Sec. 3.1, the hard-wall 
model predicts A^ ~ 2n -|- L, in contrast to the usual Regge behavior Ai^ r^ n + L found 
in experiment [50]. The radial modes are also not well described in the truncated-space 
model. For example, the first AdS radial state has a mass 1.85 GeV, which is thus difficult 
to identify with the Roper A^(1440) resonance. This problem is not present in the soft wall 
model for baryons discussed below. 

5.2 A soft- wall model for baryons 

For fermion fields in AdS one cannot break conformality with the introduction of a dilaton 
in the action since it can be rotated away leaving the action conformally invariant. ^^ As a 
result, one must introduce an effective confining potential V{z) in the action of a Dirac field 
propagating in AdS^+i space to break the conformal invariance of the theory and generate 
a baryon spectrum 

Sf = j d'^xdz^ (U'e'iT^DM^ -'-{DM^)efTH - ^i^^ -V{z)^A . (96) 

The variation of the action (96) leads to the Dirac equation in AdS 

I (^r/^'^FM^jv + ^F, j - /ii? - RV{z) * = 0. (97) 

As in the case for the hard wall model described in the previous section, the corresponding 
light-front wave equation in physical space-time follows from identifying the transverse LF 
coordinate C with the AdS holographic variable 2;, 2 — >■ (^, and the substitution (88) in (97). 
For d = 4 we find the matrix eigenvalue equation in the 2x2 spinor component representation 

-^^+ + ^^+ + t/(C)V^+ = -MV'-, 
dC C 

--^^_ + ^V^„ + t/(C)^_ = M^+, (98) 

dC C 

where U{Q = f V{Q is the effective confining potential in the light-front Dirac equation. 



^^This remarkable property was first pointed out in Ref. [132], and later derived independently in Ref. [133]. 
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Instead of choosing a dilaton profile to reproduce linear Regge behavior, as described in 
Sec. 3.2 for the case of mesons, we choose the confining interaction V in (96) to reproduce 
linear Regge trajectories for the baryon mass spectrum A^^. This "soft-wall" model for 
baryons in a higher dimensional AdS space, has also a LF analogue; it corresponds to a 
Dirac equation in physical space-time in presence of an effective linear confining potential 
U defined at equal LF time. For the potential U = k'^( equation (98) is equivalent to the 
system of second order equations 

"^^ ^~^''\ n\'^ + 2{u + 1)A MO = M'MC), (99) 



and 



^ ^ ^^''^^^- + kV + 2u^']^P-{C)=M'MC). (100) 



dC 4(2 

As a consequence, when one squares the Dirac Equation with U{(), one generates a Klein- 
Gordon equation with the potential k'^z'^. This is consistent with the same confining potential 
which appears in the meson equations. The LF equation Hlf'4'± = M.'^ip± has thus the two- 
component solution 

^^+(0 - C'^^'^e-'^'/'Ki^'e), ^40 - C^+^e-^^^^/2Lr^(K^C^), (101) 

with equal probability for the properly normalized components. The eigenvalues are 

A^2^ V(n + z/ + l), (102) 

identical for both plus and minus eigenf unctions. Note that, as expected, the potential k^C^ 
in the second order equation matches the soft-wall potential for mesons discussed in Sec. 
3.2. However, in contrast to the case for mesons, the dilaton modification of the action gives 
little guidance for finding an effective potential for baryons, since the dilaton can be scaled 
away by a field redefinition. Consequently the overall energy scale is left unspecified for 
the baryons [121]. The remarkable regularities observed in the nucleon spectrum and the 
analytical properties of the AdS/LF equations allows us, nonetheless, to built precise rules 
to describe the observed baryon spectrum and make predictions for, as yet undiscovered, 
new baryon excited states. 

Before computing the baryon spectrum we must fix the overall mass scale and the param- 
eter z/. Since our starting point for finding the bound state equation of motion for baryons is 
the light-front method, we shall require the mass scale to be identical for mesons and baryons 
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Figure 8: Orbital and radial baryon excitations for the positive-parity Regge trajectories for the 
A^ (left) and A (right) families for k = 0.49 - 0.51 GeV. 

while maintaining chiral symmetry for the pion [121] in the LF Hamiltonian equations. In 
practice, these constraints require a subtraction of — 4k^ from (102). ^^ 

As is the case for the truncated-space model, the value of z/ is determined by the short 
distance scaling behavior, u = L + 1. Higher-spin fermionic modes '$^-^...^j_^,^, J > 1/2, with 
all of its polarization indices along the 3 + 1 coordinates follow by shifting dimensions for 
the fields as shown for the case of mesons in Ref. [54] ^^. Therefore, as in the meson sector, 
the increase in the mass A4^ for baryonic states for increased radial and orbital quantum 
numbers is An = Ak"^, AL = 4k^ and AS = 2k^, relative to the lowest ground state, the 
proton; i.e., the slope of the spectroscopic trajectories in n and L are identical. Thus for the 
positive-parity nucleon sector 



M 



2(+) 
n,L,S 



Ak' 



^ S 3 

n + LH h - 

2 4 



(103) 



ior ^ 

2 ^^ 2- 



where the internal spin S 

The resulting predictions for the spectroscopy of positive-parity light baryons are shown 
in Fig. 8. Only confirmed PDG [49] states are shown. The Roper state A^(1440) and 



^^This subtraction to the mass scale may be understood as the displacement required to describe nucleons 
with Nc = 3 as a composite system with leading twist 3 + L; i.e., a quark-diquark bound state with a twist-2 
composite diquark rather than an elementary twist- 1 diquark. 

^•^The detailed study of higher fermionic spin wave equations in modified AdS spaces is based on our 
collaboration with Hans Gucntcr Dosch [32]. Sec also the discussion in Rcf. [33]. 
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the A^(1710) are well accounted for in this model as the first and second radial states of 
the proton. Likewise, the A(1660) corresponds to the first radial state of the A(1232) as 
shown in in Fig. 8. The model is successful in explaining the parity degeneracy observed in 
the light baryon spectrum, such as the L = 2, A^(1680) — A^(1720) degenerate pair and the 
L = 2, A(1905), A(1910), A(1920), A(1950) states which are degenerate within error bars. 
The parity degeneracy of baryons shown in Fig. 8 is also a property of the hard-wall model 
described in the previous section, but in that case the radial states are not well described [51]. 

In order to have a comprehensive description of the baryon spectrum, we need to extend 
(103) to the negative-parity baryon sector. In the case of the hard-wall model, this was 
realized by choosing the boundary conditions for the plus or minus components of the AdS 
wave function \1/^. In practice, this amounts to allowing the negative-parity spin baryons to 
have a larger spatial extent, a point also raised in [134]. In the soft- wall model there are no 
boundary conditions to set in the infrared since the wave function vanishes exponentially for 
large values of z. We note, however, that setting boundary conditions on the wave functions, 
as done in Sec. 5.1, is equivalent to choosing the branch u = fiR — | for the negative- 
parity spin-i baryons and u = fiR + | for the positive parity spin-| baryons. This gives 
a factor Ak,"^ between the lower-lying and the higher-lying nucleon trajectories as illustrated 
in Fig. 9, where we compare the lower nucleon trajectory corresponding to the J = L + S 
spin-| positive-parity nucleon family with the upper nucleon trajectory corresponding to the 
J = L + S — 1 spin- 1 negative-parity nucleons. As is clearly shown in the figure, the gap is 
precisely the factor 4k^. 

If we apply the same spin-change rule previously discussed for the positive-parity nucle- 
ons, we would expect that the trajectory for the family of spin-| negative-parity nucleons 
is lower by the factor 2k^ compared to the spin-| minus-parity nucleons according to the 
spin-change rule previously discussed. Thus the formula for the negative-parity baryons 

-^'ti = 4«^' (^ + ^ + f + i) ' (104) 

where S* = | or |. It is important to recall that our formulas for the baryon spectrum are 
the result of an analytic inference, rather than formally derived. 

The full baryon orbital excitation spectrum listed in Table 2 for n = is shown in Fig. 
10. We note that A^^^ c_3 = -^ni s-- ^^^ consequently the positive and negative-parity A 
states lie in the same trajectory, consistent with the experimental results. Only the confirmed 
PDG [49] states listed in Table 2 are shown. Our results for the A states agree with those 
of Ref. [59]. "Chiral partners" as the A^(1535) and the A^(940) with different orbital angular 
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Figure 9: Spectrum gap between the negative-parity spin-| nucleons and the spin-i positive-parity 



nucleon famihes for k = 0.49 GeV. 



momentum are non-degenerate from the onset. Using (103) and (104) we find the relation 

A^Af(1535) 



2' 



(105) 



A^Af(940) 

which is consistent with experiment to a good accuracy. One can in fact also build the 
entire negative-parity excitation spectrum starting from the proton partner, the J = 1/2 
negative-parity nucleon state A^(1535), using the same rules e.g., an increase in mass Ai^ 
of 4/t;^ for a unit change in the radial quantum number, 4/€^ for a change in one unit in 
the orbital quantum number and 2k'^ for a change of one unit of spin relative to the lowest 
negative-parity state, the A^(1535). 

With the exception of the A(1930) state (which is not included in Table 2), all the 
confirmed baryon excitations are well described by formulas (103) and and (104). If we follow 
the non-SU{6) quantum number assignment for the A(1930) given in Ref. [130], namely 
S = 3/2, L = 1, n = 1 we find from (104) the value A^a{1930) = 4^ ^ 2 GeV, consistent with 
the experimental result 1.96 GeV [49]. Expected results from new experiments are important 
to find out if new baryonic excitations follow the simple pattern described by Eqs. (103) and 
(104). 

An important feature of light-front holography is that it predicts a similar multiplicity 
of states for mesons and baryons, consistent with what is observed experimentally [50]. This 
remarkable property could have a simple explanation in the cluster decomposition of the 
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Figure 10: Baryon orbital trajectories for the A^ (left) and A families (right) for n = and 
K = 0.49 — 0.51 GeV. The lower and upper nucleon trajectories (left) correspond respectively to 
the the spin-2, positive-parity, and to the spin-|, negative-parity, families. The middle dotted 
trajectory (left) corresponds to spin-2 negative-parity nucleons. Plus and minus-parity states for 
the A states (right) are in the same Regge trajectory. 

holographic variable, which labels a system of partons as an active quark plus a system of 
n — 1 spectators. From this perspective, a baryon with n = 3 looks in light-front holography 
as a quark-scalar-diquark system. It is also interesting to notice that in the hard wall model 
the proton mass is entirely due to the kinetic energy of the light quarks, whereas in the 
soft-wall model described here, half of the invariant mass squared A^^ of the proton is due 
to the kinetic energy of the partons and half is due to the confinement potential. 

6 Nucleon form factors 

Proton and neutron electromagnetic form factors are among the most basic observables of 
the nucleon, and thus central for our understanding the nucleon's structure and dynamics. 
In general two form factors are required to describe the elastic scattering of electrons by 
spin-^ nucleons, the Dirac and Pauli form factors, Fi and F2 



[pv^mp) = u{p') 



rF,{q' 



2M 



F2{q' 



u{P), 



(106) 



where q = P' — P. In the light-front formalism one can identify the Dirac and Pauli form 
factors from the LF spin-conserving and spin-flip current matrix elements of the J~^ cur- 
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rent [135]. 

On the higher dimensional gravity side the spin-non-flip amplitude for the EM transition 
corresponds to the non-local coupling of an external EM field A^'^{x, z) propagating in AdS 
with a fermionic mode \l/p(a;, 2), given by the left-hand side of the equation below 

I (fx dz ^^! P^x, z) eiiV aA^ {x, z)^! p{x, z) 

~ (27r)^<5^ {P'-P- q) e,u{P')rF,{q')u{P), (107) 

where efj = (-j) 6fj is the vielbein with curved space indices M,N = 1, ■ ■ -5 and tangent 
indices A,B = 1,- ■■ ,5. The expression on the right-hand side represents the Dirac EM 
form factor in physical space-time. It is the EM matrix element (106) of the local quark 
current J^ = Cgq'j'^q with local coupling to the constituents. In this case one can also show 
that a precise mapping of the J"*" elements can be carried out at fixed LP time, providing 
an exact correspondence between the holographic variable z and the LP impact variable ( 
in ordinary space-time with the result [31] 

G±m=9±R' j -,V{Q\z)ml{z), (108) 

for the components \1/+ and \l/_ with angular momentum L^ = and L^ = +1 respectively. 
The effective charges g+ and g^ are determined from the spin-flavor structure of the theory. 
A precise mapping for the Pauli form factor using light-front holographic methods has not 
been carried out. To study the spin-flip nucleon form factor F2 using holographic methods, 
Abidin and Carlson [136] propose to introduce a non-minimal electromagnetic coupling with 
the 'anomalous' gauge invariant term 

J d^x dz^^ ei, el [Ta, Tb] F*^^v|/, (109) 

in the five- dimensional action, since the structure of (107) can only account for Fi. Although 
this is a practical avenue, the overall strength of the new term has to be fixed by the static 
quantities and thus some predictivity is lost. 

Light-front holographic QCD methods have also been used to obtain generalized parton 
distributions (GPDs) of the nucleon in the zero skewness limit in Refs. [137] and [138] for 
the soft and hard-wall models respectively. GPDs are nonperturbative, and thus holographic 
methods are well suited to explore their analytical structure. ^^ In the sections below we 
discuss elastic and transition nucleon form factors using light-front holographic ideas. ^^ ^^ 



^''^See also the discussion in Ref. [139]. 

^^A study of the EM nucleon to A transition form factors has been carried out in the framework of the 
Sakai and Sugimoto model in Ref. [140]. 

^^LF holographic methods can also be used to study the flavor separation of the elastic nucleon form 
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6.1 Computing nucleon elastic form factors in light-front holo- 
graphic QCD 

In order to compute the individual features of the proton and neutron form factors one 
needs to incorporate the spin-flavor structure of the nucleons, properties which are absent 
in models of the gauge/gravity correspondence. The spin-isospin symmetry can be readily 
included in AdS/QCD by weighting the different Fock-state components by the charges and 
spin-projections of the quark constituents; e.g., as given by the SU{6) spin-flavor symmetry. 
We label by Ng-^ and A^^^^ the probability to find the constituent g in a nucleon with spin up 
or down respectively. For the SU{6) wave function we have 

iV.t = | N^i = l, ^^^ = 1^ N,^ = l, (110) 

for the proton and 

Nut = l, ^«^ = 3' ^-^^=3' ^^^^=3' (111) 

for the neutron. The effective charges g+ and g^ in (108) are computed by the sum of the 
charges of the struck quark composed by the corresponding probability for the L^ = and 
L^ = +1 components \&+ and \l/_ respectively. We find Qp = 1, Qp = 0, g^ = —■^ and 
gf" = |. The nucleon Dirac form factors in the SU{6) limit are thus given by 

F?m = R'I^V{Q\z)^l{z), (112) 

F-m = -^-R^I^V{Q',z)[^l{z)-^l{z)], (113) 

where Ff (0) = 1 and Ff (0) = 0. 

In the soft-wall model the plus and minus components of the twist-3 nucleon wave func- 
tion are 

^+(^) = ^z'/'e-'^'/\ ^.{z) = ^,z'"e-'^'^''\ (114) 

and V{Q'^,z) is given by (60). The results for Ff'" follow from the analytic form (63) for 
any twist r. We find 

Ff(Q2) = F+(Q2), (115) 

and 

F^m = -\{F+{Q')-F^{Q')): (116) 



factors which have been determined recently up to Q^ ~ 3.4 GcV^ [141]. This wiU be described elsewhere. 
See also Ref. [142]. 
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Figure 11: Predictions for Q'^Ff (Q^) (left) and Q^Ff (Q^) (right) in the soft wah model. Data 
compilation from Diehl [143]. 



where we have, for convenience, defined the twist-2 and twist-3 form factors 

FAQ') = ^ 






and 



F-{Q' 



1 _L J^ 



Mijv ' M^,) V ^" 'mf;, 



Q2 



:ii7) 



118) 



As discussed in Sec. 4.2, the multiple pole structure in (117) and (118) is derived from the 
dressed EM current propagating in AdS space. 

The results for Q'^FfiQ'^) and Q^F^IQ^) are shown in Fig. 11. To compare with physical 
data we have shifted the poles in expression (63) to their physical values located at M^ = 
4:K?{n + 1/2) following the discussion in Sec. 4.4. The value k = 0.545 GeV is determined 
from the p mass. 

The expression for the elastic nucleon form factor Ff'" follows from (106) and (109). 



lP,r^//-)2^ 



FriQ 



dz 



^!+{z)V{Q\z)^_{z). 



(119) 



Using the twist-3 and twist-4 AdS soft-wall wavefunctions \1/_|_ and \l/_ (114) we find 

Frm = Xv.nF.{Q'), (120) 
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Figure 12: Predictions for Fj(Q^) (left) and F^(Q^) (riglit) in the soft wall model. Data compi- 
lation from Diehl [143]. 



where the amplitude (119) has been normalized to the static quantities Xp ^-^d Xn and F^{Q^) 
is given by (118). The experimental values Xp = 1-793 and Xn = —1-913 are consistent with 
the SU{6) prediction [144] ^p/fiN = —3/2. In fact (/ip//i Ar)oxp = —1-46 where /ip = 1 + Xp 
and fij\[ = Xn- The results for F|'((5^) and F2{Q^) for k = 0.545 GeV are shown in Fig. 12. 
We compute the charge and magnetic root-mean-square (rms) radius from the usual 
electric and magnetic nucleon form factors 



Q 



GE{q') = F,{q') + -^Mr) 



and 



Using the definition 



GM{q' 



{r') 



AM' 
- F,{q') + F,{q'] 
6 dF{Q^) 



F(0) rfQ2 



(121) 



(122) 



(123) 



we find the values ^/{rE)p = 0.802 fm, ^/(rl^ = 0.758 fm, (r|)„ = -0.10 fm^ and 
\/ {r\j)n = 0.768 fm, compared with the experimental values \/{rE)p = (0.877 ± 0.007) fm, 
^/{ri)p = (0.777 ± 0.016) fm, (r|)„ = (-0.1161 ± 0.0022) fm^ and v^K/)™ = (0.862 ± 
0.009) fm from electron-proton scattering experiments [49]. ^^ The muonic hydrogen mea- 
surement gives \/{rE)p = 0.84184(67) fm from Lamb-shift measurements [145]. ^^ 
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The neutron charge radius is defined by (r|;) 



_ f. dGE(Qn 



32=0 



^^Other soft and hard- wall model predictions of the nucleon rms radius are given in Refs. [136, 137, 138]. 
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Chiral effective theory predicts that the slopes are singular for zero pion mass. For 
example, the slope of the Pauli form factor of the proton at g^ = computed by Beg 
and Zepeda diverges as 1/m^ [146]. This comes from the simple triangle diagram 7* — )■ 
TT+TT^ — !■ pp. One can also argue from dispersion theory that the singular behavior of the 
form factors as a function of the pion mass comes from the two-pion cut. Lattice theory 
computations of nucleon form factors require in fact the strong dependence at small pion mass 
to extrapolate the predictions to the physical pion mass [147]. The two-pion calculation [146] 
is a Born computation which probably does not exhibit vector dominance. To make a reliable 
computation in the hadronic basis of intermediate states one evidently has to include an 
infinite number of states. On the other hand, chiral divergences do not appear in AdS/QCD 
when we use the dressed current since, as shown is Sec. 4.2, the holographic analysis with 
a dressed EM current in AdS generates instead a nonperturbative multi-vector meson pole 
struct ure.^^ 



6.2 Computing nucleon transition form factors in light-front holo- 
graphic QCD 

As an illustrative example we consider in this section the form factor for the 7*p — ?■ 
A^(1440)Pii transition measured recently at JLab. We shall weight the different Fock-state 
components by the charges and spin-projections of the quark constituents using the SU{6) 
spin-flavor symmetry as in the previous section. The expression for the spin non-flip proton 
form factors for the transition n, L — )■ n'L is [31] 

F?n,L^n',L{Q') = R' J ^ ^^ \z)V {Q', z)^-^ \z) , (124) 

where we have factored out the plane wave dependence of the AdS fields 



n^+^ I 2n\ 



The orthonormality of the Laguerre polynomials in (125) implies that the nucleon form 
factor aX Q^ = is one if n = n' and zero otherwise. Using the integral representation of 
the bulk-to-boundary propagator V{Q'^,z) given by (62) we find the twist-3 spin non-flip 



^^In the limit of a free propagating current in AdS, we obtain logarithmic divergent results: {r'i)Fi 
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Figure 13: Proton transition form factor -^f jv->Af' 
JLAB [148]. 



(Q^) to the first radial excited state. Data from 



transition form factor 



^IN^N'iQ 



V2 



Ml 
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1 + 



Q2 



(126) 



The result (126), compared with available data in Fig. 13, correspond to the valence 
approximation. The transition form factor (126) is expressed in terms of the mass of the 
p vector meson and its first two radial excited states, with no additional parameters. The 
results in Fig. 13 are in good agreement with experimental data. The transition to the 
A^(1440)Pii state corresponds to the first radial excitation of the three-quark ground state 
of the nucleon. In fact, the Roper resonance iV(1440)Pii and the A^(1710)Pii are well 
accounted in the light-front holographic framework as the first and second radial states of 
the nucleon family as shown in Sec. 5.2 (See Fig. 8). It is certainly worth to extend the simple 
computations described here and perform a systematic study of the different transition form 
factors measured at JLab. This study will help to discriminate among models and compare 
with the new results expected from the JLab 12 GeV Upgrade Project, in particular at 
photon virtualities Q^ > 5 GeV^, which correspond to the experimental coverage of the 
CLAS12 detector at JLab [149]. 
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7 Higher Fock components in light-front holographic 
QCD 

The light-front Hamiltonian eigenvalue equation (7) is a matrix in Fock space which rep- 
resents an infinite number of coupled integral equations for the Fock components ipn = (^l"^)- 
The resulting potential in quantum field theory can be considered as an instantaneous four- 
point effective interaction in LF time, similar to the instantaneous gluon exchange in the 
light-cone gauge A'^ = 0, which leads to qq — )■ qq, qq -^ qq, q -^ qqq and q -^ qqq as 
in QCD(1+1). Higher Fock states can have any number of extra qq pairs, but surpris- 
ingly no dynamical gluons. Thus in holographic QCD, gluons are absent in the confinement 
potential. ^° This unusual property of AdS/QCD may explain the dominance of quark inter- 
change [152] over quark annihilation or gluon exchange contributions in large angle elastic 
scattering [153]. ^^ 

In order to illustrate the relevance of higher Fock states and the absence of dynamical 
gluons at the hadronic scale, we will discuss a simple semi-phenomenological model of the 
elastic form factor of the pion where we include the first two components in a Fock expansion 
of the pion wave function |7r) = 4'qq/TT\qq)T=2 + '^qqqq\'i'i(l'i)r=4: + ■ ■ ■ , where the J^^ = 0~^ 
twist-two and twist-4 states \qq) and \qqqq) are created by the interpolating operators q^'^'^^q 
and q'j^'j^qqq respectively. 

Since the charge form factor is a diagonal operator, the final expression for the form factor 
corresponding to the truncation up to twist four is the sum of two terms, a monopole and 
a three-pole term. In the strongly coupled semiclassical gauge/gravity limit hadrons have 
zero widths and are stable. One can nonetheless modify the formula (63) by introducing 
a finite width: q^ ^ q^ + ^/2iMT. We choose the values F^ = 140 MeV, F^/ = 360 MeV 
and Tpii = 120 MeV. The results for the pion form factor with twist two and four Fock 
states are shown in Fig. 14. The results correspond to Pqqqq = 13 %, the admixture of the 
\qqqq) state. The value of Pqqqq (and the widths) are input in the model. The value of k 
is determined from the p mass and the masses of the radial excitations follow from setting 
the poles at their physical locations, A^^ — )■ An^iji + 1/2), as discussed in Sec. 4.4. The 
time-like structure of the pion form factor displays a rich pole structure with constructive 



•^^This result is consistent with the flux-tube interpretation of QCD [150] where soft gluons interact so 
strongly that they are sublimated into a color confinement potential for quarks. The absence of constituent 
glue in hadronic physics has been invoked also in Ref. [151], where the role of the confining potential is 
attributed to an instanton induced interaction. 

■^^In Rcf. [154] we discuss a number of experimental results in hadron physics which support this picture. 
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Figure 14: Structure of the space-like (left) and time-like (right) pion form factor in light-front 
holography for a truncation of the pion wave function up to twist four. Triangles are the data 
compilation from Baldini et al. [77], squares are JLAB data [78]. 

and destructive interfering phases; this is incompatible with the admixture of the twist-three 
state \qqg) containing a dynamical gluon since the interference in this case is opposite in 
sign. 



8 Conclusions 

As we have shown, the exact light-front Hamiltonian Hlp\iI)) = A^^|^) for QCD can be 
systematically reduced to a relativistic frame-independent semiclassical wave equation [12] 



d^ 



AL' 



de 



<' 



t/(C)U(C) = -MV(C), 



(127) 



for the valence Fock state of mesons. The unmodified AdS equations correspond to the kinetic 
energy terms of the massless constituent quarks with relative orbital angular momentum L = 
L^. The effective potential U{() corresponds to the color-confining potential and follows from 
the truncation of AdS space, in a modified effective AdS action, and light-front holography. 
The variable ( is the invariant separation of the constituents. This frame-independent light- 
front wave equation is comparable in simplicity to Schrodinger theory in atomic physics 
which is formulated at equal instant time. We have also derived an analogous light-front 
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Dirac equation for holographic QCD which describes hght-quark baryons with finite color 

Remarkably, these light-front equations are equivalent to the equations of motion in 
a higher dimensional warped space asymptotic to AdS space. The mapping of the gravity 
theory to the boundary quantum field theory, quantized at fixed light-front time, thus gives a 
precise relation between holographic wave functions and the light-front wave functions which 
describe the internal structure of the hadrons and their electromagnetic couplings. This 
mapping provides the basis for a profound connection between physical QCD quantized 
in the light-front and the physics of hadronic modes in a higher dimensional AdS space. 
However, the derivation of the effective color-confining potential ^(C) directly from QCD, 
remains an open question. 

Despite some limitations of AdS/QCD [155], the light- front holographic approach to 
the gauge/gravity duality, Light-Front Holography, has already provided significant physical 
insight into the strongly-coupled nature and internal structure of hadrons; in fact, it is one of 
the few tools available. As we have seen, the resulting model provides a simple and successful 
framework for describing nonperturbative hadron dynamics: the systematics of the excitation 
spectrum of hadrons: the mass eigenspectrum, observed multiplicities and degeneracies. It 
also provides powerful new analytical tools for computing hadronic transition amplitudes, 
incorporating the scaling behavior and the transition from the hard-scattering perturbative 
domain, where quark and gluons are the relevant degrees of freedom, to the long range 
confining hadronic region. 

The dressed current in AdS includes the nonperturbative pole structure. Consequenly, 
the approach incorporates both the long-range confining hadronic domain and the constituent 
conformal short- distance quark particle limit in a single framework. The results display a 
simple analytical structure which allows us to explore dynamical properties in Minkowski 
space-time; in many cases these studies are not amenable to Euclidean lattice gauge theory 
computations. In particular, the excitation dynamics of nucleon resonances encoded in the 
nucleon transition form factors can provide fundamental insight into the strong-coupling 
dynamics of QCD. New theoretical tools are thus of primary interest for the interpretation 
of the results expected at the new mass scale and kinematic regions accessible to the JLab 
12 GeV Upgrade Project. 

The semiclassical approximation to light-front QCD described in this article is expected 
to break down at short distances where gluons become dynamical degrees of freedom and 
hard gluon exchange and quantum corrections become important. One can systematically 
improve the semiclassical approximation, for example, by introducing nonzero quark masses 
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and short-range Coulomb-like gluonic corrections, thus extending the predictions of the 
model to the dynamics and spectra of heavy and heavy-light quark systems. The model 
can also be improved by applying Lippmann-Schwinger methods to systematically improve 
the light-front Hamiltonian of the semiclassical holographic approximation. One can also 
use the holographic LFWFs as basis functions for diagonalizing the full light-front QCD 
Hamiltonian[156] as well as the input boundary functions to study the evolution of structure 
functions and distribution amplitudes at a low energy scale. 
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Appendices 



A AdS boundary conditions and interpolating opera- 
tors 

The formal statement of the duahty between a gravity theory on {d + l)-dimensional 
Anti-de Sitter AdSd+i space and the strong couphng hmit of a conformal field theory (CFT) 
on the d-dimensional asymptotic boundary of AdSd+i at z = is expressed in terms of the 
d+ 1 partition function for a field $(a;, z) propagating in the bulk 

Zgra.m = e''^"^"^ = /"l?[<l>]e^^[*l, (128) 

where Se/f is the effective action of the AdSd+i theory, and the d- dimensional generating 
functional of correlation functions of the conformal field theory in presence of an external 
source ^o{x'^) 

Zcft[%] = e^^^^^[*«l = /exp (i f d'^x $o(x)C(x) J \ . (129) 

The functional Wcft is the generator of connected Green's functions of the boundary theory 
and (9 is a QCD local interpolating operator. 

According to the AdS/CFT correspondence, to every operator in the conformal field 
theory there corresponds an AdS field. We use the isometrics of AdS space to map the 
scaling dimensions of the local interpolating operators defined at the AdS boundary into 
the modes propagating inside AdS space. The precise relation of the gravity theory on AdS 
space to the conformal field theory at its boundary is [7] 

Zgrav[^{x,z)\^^^ = %{x)\ =ZcFt[%], (130) 

where the partition function (128) on AdSd+i is integrated over all possible configurations 
$ in the bulk which approach its boundary value $o- If we neglect the contributions from 
quantum fluctuations to the gravity partition function, then the generator Wcft of connected 
Green's functions of the four-dimensional gauge theory (129) is precisely equal to the classical 
(on-shell) gravity action (128) 

Wcft [M = S,ff[^{x,z)\^^^ = $o(a:)]„^,_,,,u, (131) 
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evaluated in terms of the classical solution to the bulk equation of motion. This defines the 
semiclassical approximation to the conformal field theory. In the bottom-up phenomeno- 
logical approach, the effective action in the bulk is usually modified for large values of z to 
incorporate confinement and is truncated at the quadratic level. 
In the limit 2; — ?■ 0, the independent solutions behave as 

<l>(x,2) ^2^$+(x)+z'^~^<l>_(x), (132) 

where r is the scaling dimension. The non-normalizable solution $„ has the leading boundary 
behavior and is the boundary value of the bulk field $ which couples to a QCD gauge 
invariant operator O in the 2; — )■ asymptotic boundary, thus $_ = $0- The normalizable 
solution $_|_ is the response function and corresponds to the physical states [157]. The 
interpolating operators O of the boundary conformal theory are constructed from local gauge- 
invariant products of quark and gluon fields and their covariant derivatives, taken at the same 
point in four- dimensional space-time in the x^ — > limit. According to (129) the scaling 
dimensions of O are matched to the conformal scaling behavior of the AdS fields in the limit 
2; — ;■ and are thus encoded into the propagation of the modes inside AdS space. 

Integrating by parts, and using the equation of motion for the field in AdS, the bulk 
contribution to the action vanishes, and one is left with a non-vanishing surface term in the 
ultraviolet boundary 



1 



S = R"-^ lim / (fx -T— $a^$, (133) 



2-s>0 ' ^'^-l 



which can be identified with the boundary QFT functional Wcft- Substituting the leading 
dependence (132) of $ near z = in the ultraviolet surface action (133) and using the 
functional relation 

one finds that ^+{x) is related to the expectation values of O in the presence of the source 

$0 [157] 

(0|O(x)|0)^„~<l>+(x). (135) 

The exact relation depends on the normalization of the fields chosen [158]. The field $+ 
thus acts as a classical field, and it is the boundary limit of the normalizable string solution 
which propagates in the bulk. 
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